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ABSTRACT 

We report the results of a 15 GHz (2 cm) multi-epoch VLBA program, begun in 1994 to study the outflow 
in radio jets ejected from quasars and active galaxies. The observed flow of 208 distinct features measured 
in 1 10 quasars, active galaxies, and BL Lac objects shows highly collimated relativistic motion with apparent 
transverse velocities typically between zero and about 15c, with a tail extending up to about 34c. Within in- 
dividual jets, different features appear to move with a similar characteristic velocity which may represent an 
underlying continuous jet flow, but we also see some stationary and even apparently inward moving features 
which co-exist with the main features. Comparison of our observations with published data at other wave- 
lengths suggests that there is a systematic decrease in apparent velocity with increasing wavelength, probably 
because the observations at different wavelengths sample different parts of the jet structure. 

The observed distribution of linear velocities is not consistent with any simple ballistic model. Either there 
is a rather broad range of Lorentz factors, a significant difference between the velocity of the bulk relativistic 
flow and the pattern speed of underlying shocks, or a combination of these options. Assuming a ballistic flow, 
comparison of observed apparent velocities and Doppler factors computed from the time scale of flux density 
variations is consistent with a steep power law distribution of intrinsic Lorentz factors, an isotropic distribution 
of orientations of the parent population, and intrinsic brightness temperatures about an order of magnitude 
below the canonical inverse Compton limit. It appears that the parent population of radio jets is not dominated 
by highly relativistic flows, and contrary to the assumption of simple unified models, not all sources have 
intrinsic speeds close to c. 

Usually, the observed jet flow is in the general direction of an established jet. However, many jets show 
significant bends and twists, where the observed motions are non-radial, but are alingned with the local jet 
direction suggesting that the jet flow occurs along pre-existing bent channels. In a few cases we have observed 
a clear change in the direction of a feature as it flows along the jet. Radio jets which are also strong gamma-ray 
sources detected by EGRET appear to have significantly faster speeds than the non EGRET sources, consistent 
with the idea that gamma ray sources have larger Doppler factors than non gamma-ray sources. Sources at high 
redshift have systematically lower angular speeds than low redshift jets, consistent with standard cosmologies. 
Subject headings: galaxies: active — quasars: general — galaxies: jets — radio continuum: galaxies 
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1. INTRODUCTION 

The discovery in the mid 1960s of rapid variability in 
extragalactic radio sources iSholom itskii 1965; Dent 1965 
iPaulinv-Toth & KellermanDlll966ir appeared difficult to ex- 
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plain in terms of conventional synchrotron radiation theory. 
These early observations showed changes by as much as 25 
percent over a few weeks as well as significant day-to-day 
variations (Paulinv-Toth & Kellermann 1966). It was quickly 
realized (e.g., Hovle 1966), that the observed rapid variabil- 
ity implied such small linear dimensions, that the relativis- 
tic electron population would be rapidly extinguished by in- 
verse Compton scattering. Later. iKellermann & Paulinv-TothI 
0969) put these arguments on a quantitative observational ba- 
sis showing that due to inverse Compton cooling, the max- 
imum sustainable brightness temperature for an incoherent 
source of el ectron synchro tron radiation is less th an 10^^ K. 

However- IWoltieiU 1 966h and lReesUl 96611 1967h pointed out 
that if there is relativistic bulk motion of the emitting mate- 
rial, since the radiation is then beamed along the direction of 
motion, the apparent luminosity in that direction is enhanced, 
while at the same time the cross-section for inverse Comp- 
ton scattering is reduced. Additionally, since the radiating 
source nearly catches up with its radiation, for a favorably po- 
sitioned observer, the apparent transverse motion can exceed 
the speed of light. Early VLBI observations (Whitnev et al. 
[T97I; Cohen et al. 1971) showed evidence for the predicted 
high-velocity outflow; however, the arguments were indirect, 
and one had to have faith in the interpretation of the very 
limited radio interferometric data which did not fully sam- 
ple the source structure (e.g.. Dent 1972). Subsequent higher 
quality VLBI observations confirmed the existence of super- 
luminal motion in the well-coUimated radio jets, found in the 
nuclei of many qu asars as well as in nearby active galaxies 
llCohenetal.11977'). 

VLBI observations provide a direct method to investigate 
aspects of the formation, acceleration, and propagation of ex- 
tragalactic radio jets. Early studies discussing apparent jet 
speeds were statistically unreliable because they were largely 
based on the analysis of published observations of only a few 
tens of sources made at different times, often at only two 
or three epochs, by different groups, using different antenna 
array configurations, and different observing/data reduction 
strategies. The dynamic range of the images was often inad- 
equate to identify individual features, especially for the more 
complex sources, and spacings in time between successive ob- 
servations were frequently too long to uniquely identify and 
track moving jet features from epoch to epoch. Moreover, the 
source selection criteria for many previously published stud- 
ies of superluminal motion were not well defined. A system- 
atic study by Vermeulen ( 1995) of 81 sources from the flux 
density limited 6 cm Caltech-Jodrell Bank Survey indicated 
smaller typical speeds than in earlier reports. This suggested 
that earlier studies preferentially tended to observe and/or re- 
port only sources in which rapid motion was detected or sus- 
pected. 

In 1994 we began a systematic 15 GHz VLBA survey of 
relativistic outflows in a sample of over one hundred quasar 
and active galaxy radio jets. Our motivation was to study 
the distributions of velocities, bending, pattern motions, ac- 
celerations, and other complexities of the jet kinematics that 
may exist, as well as changes in the strength and morphology 
of features as they propagate along the jet. These observa- 
tions have provided the homogeneity, resolution and dynamic 
range to reliably distinguish and identify individual compo- 
nents between different epochs. The central VLBA antennas 
provide short interferometric spacings that have allowed us to 
a) track features moving for some distance down the jet where 
they become more diffuse, and b) observe the continuous jet 



rather than just the bright features which are often referred to 
as "components" or "blobs". 

We chose 15 GHz as a compromise between achieving the 
best angular resolution, and the better sensitivity and immu- 
nity from weather conditions found at lower frequencies. At 
this observing frequency the resolution of the VLBA is ap- 
proximately 0.5 mas in right ascension and between 0.5 and 
1 mas in declination. We therefore have sufficient angular res- 
olution in many cases to resolve the 2-dimensional jet trans- 
verse to its flow, which allows for tests of some theoretical 
models. Also, individual features can often be recognized at 
significantly smaller separations from the origin. On the other 
hand, our observations are less sensitive to structure and mo- 
tions of the more diffuse features located far out along the jet 
or within 0.5 mas of the core. 

Our 15 GHz data represent a significant improvement over 
previous studies of superluminal motion in AGN. First, our 
sample of sources is large, and the sample membership was 
based on criteria other than observed superluminal motion. 
Second, the speeds are better determined because a) our data 
are well-sampled and span a longer time period, and b) our 
higher image resolution res ults in less blending of f eatures. 

In Paper I of this series llKellermann et alJll998 a). we de- 
scribed the parsec-scale structure of 132 of the strongest 
known radio jets based on a single epoch of observation. In 
Paper II (Zensus et al. 2002) we discussed the structure of an 
additional 39 sources. Contour maps of all of our multi-epoch 
observations are available on our web site * along with kine- 
matic and other data on each source. In this paper we discuss 
our multi-epoch observations made through 15 March 2001 
and the kinematics derived from these observations. 

Some prelimi nary results of our program have already 
been pubHshed (IKellermann et al."1999'. '2000^■ 'Kell ermannl 
2002; Rosetal. 2002; Cohen et al. 2003; Kellerman rTet alJ 
I2003t iKovalevi i2003c iZensus et ah i2003i) . We have also 
discussed the kinematics of seve ral individual sources 
including PKS I345-H12 5 jLister e t al. 2003), NGC 1052 
JVermeulen et alJ2003a^ and 3C 279 ( H oman et al._2003j) . 

In §|2]we describe our sample and source selection criteria, 
and in §| 3|we discuss the details of our observational program. 
InSO we summarize the predictions of relativistic beaming 
models while in §14. 21 we present the observed jet kinematics. 
In § l4.3l and § l4.4l we discuss the implications ofour observa- 
tions for the nature of the relativistic flow. In § 14.51 we discuss 
the angular velocity - redshift relation, and in §|5]we summa- 
rize our results and describe our on-going observing program. 

Throughout this paper we use the following cosmological 
parameters: Hq = 70 kms~' Mpc~', fi^ = 0.3, and fl\ = 0.7. 

2. SAMPLE DEFINITION 

2.1. The Full Sample 

Our original source sample (see Paper I) was based on 
the Kiihr I Jy catalog (Kiihr et al. 1981) as supplemented by 
Stickel, Meisenheimer, & Kiihr ( 1994). Our goal was to in- 
clude all known sources in the Stickel et al. catalog that have 
a flat spectral index (a > -0.5 for 5^ ^ i^") anywhere above 
500 MHz, and a total flux density at 15 GHz, observed at least 
at one epoch, greater than 1 .5 Jy for sources north of the ce- 
lestial equator, or greater than 2 Jy for sources between decli- 
nations -20° and the equator Since the Stickel et al. catalog 
is complete only at 5 GHz, we used other measurements at 

http://www.nrao.edu/2cmsurvey/ 
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15 GHz or extrapolations from lower frequencies to augment 
the Kiihr sample. 

While we attempted to be complete for flat-spectrum, core- 
dominated sources which met our selection criteria, we also 
included a number of other sources of special interest as fol- 
lows: 

a) Six compact steep spectrum (CSS) sources as indicated 
in Table [2 These sources generally have steep radio spectra 
and angular sizes smaller than ^ 1 " on the sky. 

b) Four lobe-dominated sources with core components 
that would satisfy our criteria if they did not also have 
strong extended steep spectrum structure. Two other sources, 
NGC 1052 (0238-084) and 3C 120 (0430H-052), also have 
prominent double lobe structure, but since they are core- 
dominated on arcsecond scales at 15 GHz, they met our spec- 
tral index selection criterion. 

c) Ten Gigahertz-peaked spectru m (GPS) sources, s ome of 
which have two-sided jet structure (iLister et al.ll2002t M. L. 
Lister et al., in preparation). However, many of the GPS 
sources we observed turned out to be merely flat-spectrum 
sources whose spectrum was temporarily dominated by a 
bright, synchrotron self-absorbed component in the jet. Our 
"peaked" classification in TableQlis given only to sources that 
to our know ledge have always met the GPS s pectral crite- 
ria given by de Vries. Barthel . & O'Deal Jl997h . For a few 
sources, our classification differs from that previously pub- 
lished (e.g., O'Dea 1998) ; however we believe that our clas- 
sification is more r obust ( :Kovalevll2004l and Y. Y. Kovalev et 
al., in preparation). 

Table[nsummarizes the properties of each source discussed 
in this paper All of these sources were introduced in Pa- 
pers I and II; however, 37 sources included in Papers I and 
II are not discussed now because they were observed only 
once or twice. These have subsequently been observed and 
will be discussed in a separate paper (E. Ros et al., in prepa- 
ration). One source listed in Paper I is gravitationally lensed 
(0218H-357), so we have chosen not to include it in our sta- 
tistical discussion. Eight other sources had no detectable jet 
above our sensitivity limit. They are labeled as "naked cores" 
in Table[2and we report no motions for these sources. Three 
sources, 0316H-162, 1328H-254, and 1328H-307 reported in Pa- 
pers I or II, are not discussed here as they have complex struc- 
ture that we could not adequately image. 

Columns 1 and 2 of Table ^ give the lAU source desig- 
nation and where appropriate a commonly used alias. Col- 
umn 3 indicates whether or not the source is a member of the 
r epres entative flux density limited sample that we describe in 
§ | 2.2I The optical classification and redshift as given mainly 
by Veron-Cettv & Veron (2001j) are given in columns 4 and 
5, respectively. In column 6 we give a radio spectral clas- 
sification for each source based on RATAN-600 telescope 
observations of broad-band instantaneous spectra from 1 to 
22 GHz as described by Kovalev et al. ( 1999) . These spec- 
tra are available on our web site. For a few sources not ob- 
served at the RATAN-600 telescope, we used published (non- 
simultaneous) radio flux densities taken from the literature. 
We consider a radio spectrum to be "flat" if any portion of its 
spectrum in the range 0.6 to 22 GHz has a spectral index flat- 
ter than -0.5 and steep if the radio spectral index is steeper 
than -0.5 over this entire region. Column 7 shows the parsec 
scale radio morphology taken from Papers I and II; in column 
8 the largest total flux density seen on any of our VLBA im- 
ages, and in column 9 an indication of whether or not the radio 
source is associated with a gamma-ray detection by EGRET. 



2.2. The Representative Flux Density Limited Subsample 

The full sample of sources described in the previous section 
is useful for investigating jet kinematics in a cross-section of 
known AGN classes. However, in order to compare with the 
theoretical predictions of relativistic beaming models, a well- 
defined sample selected on the basis of beamed (not total) flux 
density is needed. Past surveys (e.g., Tavlor et al. 1996; Liste3 
2001) have attempted to accomplish this by means of a spec- 
tral flatness criterion. However, we found that this method 
eliminates some lobe-dominated active galaxies such as those 
described in § 12.11 Also, we found that the extrapolated 
15 GHz flux density based on non-simultaneous lower fre- 
quency measurements was often grossly in error, because of 
spectral curvature or variability. 

We have therefore assembled a flux density-limited sub- 
sample from the full 15 GHz VLBA survey by using our mea- 
sured VLBA flux densities as the main selection criterion. All 
sources that had a total CLEAN VLBA flux density exceed- 
ing 1.5 Jy (2 Jy for southern sources) at any epoch since 1994 
are included in this sub-sample. We excluded any sources that 
were observed on at least four occasions and never exceeded 
this limit. For survey sources with fewer than four VLBA 
epochs, we estimated the VLBA flux density at various epochs 
during this period using the source compactness and data from 
the flux densi ty monitoring pro grams at the RATAN-600 ra- 
dio telescope (i Kovalev et alJl99 9) or the Uni versity of Michi- 
gan Radio Astronomy Observatory^ ( Afler et al.ll992ll2003h . 

It is important to note that this sub-sample, although 
flux density-limited, is not complete; that is, there are ad- 
ditional compact sources that fulfill our selection criteria. 
This is partly due to the lack of a complete all-sky sur- 
vey at 15 GHz, and also to the variable nature of AGN. 
We have therefore identified a list of candidate sources 
from other recent high-frequency radio surveys, made since 
our original list w as compiled in 19 94. These include 
the WMAP survey (^Bennett et al."2003l), the VLBA calibra- 
tor survey (Beasley et al. 2002), RATAN-600 observations 
('Kovalev et al. 1999), and the high-frequency peaker survey 
(Dallacasa et al. 2000). We have made new VLBA observa- 
tions of these sources in order to assemble a complete flux- 
density limited, core-selected sample. There are 133 sources 
in this complete flux density limited sample, but so far we 
have multiple epoch observations and derived speeds for only 
71 of these sources which we define as the representative sub- 
sample. For the purposes of the statistical analysis presented 
here, we consider our present sub-sample to be representa- 
tive of a complete sample, since the general properties of 
the missing sources should not be substantially different from 
the whole. We have compared the 15 GHz luminosity dis- 
tributions of our representative sub-sample and the missing 
sources from the full sample using a Kolmogorov-Smirnov 
(K-S) test, and find no significant difference between the two 
samples. 

Although our sub-sample selection method is somewhat 
complex, it is based on the directly-measured compact flux 
density, and does not use an often unreliable spectral index 
criterion. Also, since the survey membership is not deter- 
mined from a single "snapshot" epoch, we are not excluding 
potentially interesting or highly variable sources simply be- 
cause they happened to be in a low state at the time of the 
original investigation. This increases the size of the sample 

^ http://www.astro.lsa.umich.edu/obs/radiotel/umrao.html 
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and the robustness of statistical tests on source properties. Of 
the 71 sources in our representative sub-sample, there are 53 
quasars, 12 BL Lac objects, and 6 galaxies. 

3. THE OBSERVATIONS AND ANALYSIS 

Our observations were made during 29 separate observing 
sessions between 1994 and 2001. Typically, we obtained im- 
ages of each source at three to seven epochs over this seven 
year period. Sources were usually observed at 6 to 18 month 
intervals. Those with known rapid motion were observed 
more frequently, while those with no observed jet or with only 
small observed motions were observed less often. Each ob- 
serving session lasted between 8 and 24 hours. In general 
we tried to observe only at elevations above 10° to minimize 
the effect of tropospheric absorption, phase errors, and exces- 
sive ground radiation. Except for sources at low declinations, 
where the hour angle coverage is restricted, we observed each 
source once per hour for six to eight minutes over a range of 
eight hours in hour angle. In general, all ten VLBA antennas 
were used for each observations, except when restricted by el- 
evation. We rarely used less than 8 antennas. Papers I and II, 
as well as our web site, give logs of the observations. 

Some additional observations were made in 1998 and 1999 
by L. I. Gurvits et al. (in preparation) as part of a different pro- 
gram to compare 15 GHz source structure measured with the 
VLBA to the 5 GHz structure measured using the Japanese 
HALCA space VLBI satelUte (Hirabayashi et al. 1998). We 
have used these data to supplement our own as the observa- 
tions were made using the same observing and data reduction 
procedures as used for the present program. These images are 
also available on our web site. 

Data reduction was done using a combination of AIPS and 
DIFMAP as described in Paper I. Each image was analyzed 
using the AIPS tasks JMFIT or MAXFIT to determine the 
relative positions of each definable feature at each epoch, and 
these positions were then used to calculate velocities relative 
to a presumed core component. We have assumed that the 
bright unresolved feature typically found at the end of so- 
called core-jet sources is the stationary core. Generally, the 
choice of the core is clear due to its high brightness tempera- 
ture and location. However, in a few cases, the location of 
the core is ambiguous, particularly in some sources where 
the component motions appear two-sided about a centrally- 
located core. As an aid in identifying components from one 
epoch to the next we examined each image for continuity in 
position, flux density, and structure. Usually, this procedure is 
more reliable than fitting independent models to the (m,v) data 
at each epoch, and with the possible exception of a few iso- 
lated cases, we believe that we have correctly cross-identified 
features as they evolve from epoch to epoch. However, for 
those sources where the jet is barely resolved from the VLBA 
core, we fit models to the data in the (m,v) plane in order to 
exploit the full interferometric resolution of the VLBA. 

In Figures 1 and|2]we show the angular motion in the jets 
of the 120 sources for which we have been able to determine 
motions with three or more epochs of observation. Of these 
120 sources, there are 110 sources with good quality data and 
measured redshifts for which we have been able to determine 
the linear velocity of at least one jet feature. These comprise 
13 active galaxies, 79 quasars, and 18 BL Lac Objects. In 
all, we have been able to determine reliable values for 208 
separate features found in these 1 10 jets. 

We have determined the radial angular speed, /i, of each 
definable jet feature using a linear least squares fit to the mea- 



sured component positions, relative to the presumed core. Our 
measured values of /x are given in Table |2l In Figures 1 
and|2]the slope of the line corresponds to the best fitting least 
squares fit to the speed which is shown in column 6 of Ta- 
ble 121 Figure 1 shows the plots for the one-sided jets and 
Figure 121 the two-sided ones. For well-defined components 
the formal uncertainty in the relative position is small, gen- 
erally less than 0.02 milliarcseconds (mas). Frequently, how- 
ever, the jets have a complex brightness distribution with re- 
gions of enhanced intensity that can brighten and fade with 
time. In some sources, instead of a well-collimated jet, there 
is a broad plume or more complex two dimensional structure. 
Not only is the centroid of these components poorly defined, 
but changes in the brightness distribution with time compli- 
cate the definition of their motion. Some features appear to 
move; while others appear stationary. Some may break up 
into two or more separate features, and it is often unclear how 
these moving features are related to the underlying relativis- 
tic flow. In some sources, the moving features are not well 
resolved from each other, especially when they are close to 
the core. These issues all affect our ability to make a unique 
identification of components from epoch to epoch. In a few 
cases where there are fast moving components, our observa- 
tions may be too widely spaced in time to uniquely define 
their speed. Ambiguities exist, but we believe that most of 
our component identifications are robust and that component 
speeds are reliably determined. 

We have based our error estimates of the angular velocity 
on the dispersion about the best-fit linear relation between the 
measured feature positions and time. The accuracy of each fit 
depends on the strength and size of the feature, the number of 
epochs, and the range of time covered by the observations. We 
note, however, that the measured dispersion of points about 
the best-fit line will overestimate the true scatter in component 
positions if there is acceleration such as observed for 3C 279 
(Homan et al. 2003). Additionally, the measured radial speed 
may be an underestimate of the true speed if there is a signifi- 
cant non-radial component to the motion. However, we show 
m § l4!23l that while non-radial motions are common, the vec- 
tor speeds do not differ significantly from the radial speeds, 
except in the case of 1548H-056. 

We use the following criteria to classify the quality of each 
measured velocity. 

• The component position is determined at four or more 
observations. 

• The component is a well-defined feature whose position 
can be unambiguously determined to a small fraction of 
the VLBA beamwidth. 

• The best fitting angular speed is determined to high ac- 
curacy, defined by an uncertainty < 0.08 mas per year 
or a significance^ > 5a. 

We then assign a "quality code" to each component motion as 
follows. 

• E (Excellent) denotes motions that satisfy all three of 
the above criteria. 

• G (Good) denotes motions that satisfy any two of the 
above criteria. 

* The > 5(T requirement is necessary to accommodate nearby sources with 
very high angular speed features which may have large absolute uncertainties. 
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Fig. 1a. — Plots showing the change in separation with time of features of one-sided jets for which we have measured a velocity from observations at three or 
more epochs. An asterisk denotes sources that were model fit in the {u,v) plane rather than in the image plane. Different symbols are used for each component as 
follows: B, filled triangle; C. filled five point star; D, filled circle; E, filled six point star; F, cross; G. five point star; H. square; I, triangle; J, circle with cross; K, 
circle with dot; L, circle. The solid lines denote the best least square linear fit to the data, and the slope represents the proper motion, /i, tabulated in Tablel2l 
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Fig. 2. — Same as Figure 1 but for two sided sources only. The points 
shown for 0238-084 (NGC 1052) are the brightest four components located 
on each side of the center of symmetry. No central component is visible for 
these sources, so the measured component positions are referred to a virtual 
center of symmetry as discussed bv .Ver meulen et al. 1 2003a). 



• F (Fair) denotes motions that satisfy only one of the 
above criteria. 

• P (Poor) denotes motions that do not satisfy any of the 
above criteria, or that the uncertainty in the fitted speed 
is > 0.15 mas per year (except for the > 5cr cases de- 
scribed above). 

These codes and the measured proper motions are tabulated 
m TablelH 

4. DISCUSSION 
4.1. Superluminal Motion and Relativistic Beaming 



We interpret the structural changes seen in the AGN of 
our sample in terms of the twin relativistic jet model of 
iBlandford &^eesl iITqtI . In this framework, the bulk veloc- 
ity of the relativistic flow, in units of the speed of light, c, is 
usually denoted as /3b, and 7b is the corresponding Lorentz 
factor defined by 7b = {\- (3^)~^l^. For a jet flow at an angle 9 
to the observer's line of sight, the Doppler factor, 5, is given 
by<5 = 7b'(l-/3bcos0)-i. 

Three important observational consequences of relativistic 
motion in synchrotron sources are: 

a) A Doppler frequency shift of the radiation, where the 
ratio of observed frequency, v, to the emitted frequency Ve is 
given by 

y/y. = 5. (1) 

b) A change in the observed transverse velocity due to the 
apparent time compression. The apparent transverse velocity, 
in units of c, is given by 

/3app = /3psin0/(l-/3pcos0), (2) 

where /3p = Vp/c is the pattern velocity. We make this 
distinction between pattern and bulk velocities because 
jet features can move with signi ficantly different velocity 
than the bulk jet flow (see, e.g., Lind & Blandford 19851 
jZensus, Cohen, & Unwin 1995; Vermeulen & Cohen 1994). 
This could be the case, for example, if the pattern mo- 
tion is due to the propagation of shocks. By study- 
ing 25 core-dominated qua sars selected from the literature, 
IVermeulen & Coheiil (11994 ') showed that the simplest model, 
which has only one value of 7 that is the same in all sources 
and also has /3b = /3p in all sources, is not tenable. They 
showed that either there must be a distribution of 7 among 
the sources; or that /3b (ip- Most likely, there is both some 
distribution of 7, and also a difference between 7p and 7b 

The apparent transverse velocity of an approaching com- 
ponent with a pattern Lorentz factor 7p, reaches a maximum 
apparent speed /3p7p when the component moves at an angle 
sin(0) = 1 /7p to the line of sight. For the corresponding com- 
ponent in the receding jet, the observed velocity is /3p/2. 

c) A change in the apparent flux density, S, of a moving 
component over its stationary value, ^o, by a factor 

5/5o = 5-^-", (3) 

where a is the spectral index and x has a value of 2 for 
a continuous jet or 3 for discrete components (see e.g., 
lUrrv & Padovani Fl99.5Y 

The Doppler factor is sharply peaked along the direction 
of motion, so sources with highly relativistic jets that happen 
to be pointed close to the line of sight will appear strongly 
boosted, and hence likely to be selected in a flux-limited sam- 
ple. Although the jets may be intrinsically two-sided, unless 
they are very close to the plane of the sky, they will appear 
highly asymmetric (i.e., one-sided) since the radiation from 
the receding jet is highly beamed away from the observer. 

Equation (|3j has the important consequence that the 
strongest, most compact radio sources we observe tend to 
have highly relativistic jets that are aligned close to the line 
of sight, and they will therefore likely display superluminal 
motion. The most probable angle for sources selected on the 
basis of bea med flux densitv alone is c lose to (27)"', where 
/3app - 7/2 ( IVermeulen & CohenlfT99l . The distribution of 
observed speeds and flux densities is dependent on the distri- 
butions of 7, redshift, and the intrinsic luminosity function. 

This is illustrated in Figure |3l where we show the predicted 
distribution of apparent velocities for three different models. 
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Fig. 3. — Predicted distributions of /3app when the sample has randomly 
oriented jets all with the same Lorentz factor 7 (dotted line), when the sample 
has sources drawn from such a jet population, but based on flux density which 
can be enhanced by Doppler boosting (solid line), and for a Doppler boosting 
with a range of Lorentz factors which favors low values of the Lorentz factor 
(dashed Une). 



If the effect of Doppler boosting is ignored; for example, if the 
observed motions are due to the propagation of shocks rather 
than actual bulk motion, then most jets will lie close to the 
plane of the sky and have apparent speeds, /3app close to f3p. If 
we take into account the effect of Doppler boosting, then, in a 
flux limited sample, with a single Lorentz factor, most sources 
are found to lie close to the line of sight, and have an apparent 
velocity, /J^pp ^ /3p7. In the typical situation which we con- 
sider, where /3p ~ 1, then in the absence of Doppler boosting, 
case most sources appear to have an apparent velocity close 
to c, whereas in the Doppler boosted case, most sources have 
an apparent velocity close to 7c. Note that in the absence of 
Doppler boosting, there are virtually no sources with appar- 
ent speeds less than (3c since this would require rare end-on 
orientations. The randomly oriented model also drops much 
more sharply with increasing (3 than the model with a range 
of Lorentz factors. 

Blandford & Konigj (112221) suggested that one of the ob- 
served jet components might be the stationary feature of the 
approaching jet at the point where it becomes optically thin, 
so that the flux density of this "core" component as well as the 
moving component is Doppler boosted. The observed overall 
spectrum is typically flat, which may due to the superposition 
of different parts of the jet, which each have different Doppler 
factors, causing the synchrotron cutoff to appear at differ- 
ent frequencies. Consideration of differential Doppler boost- 
ing led to the concept of unified models (e.g.. lOrr & Brown^ 
[1982; Barthel 1989; Urrv & Padovani 1995) which is usually 
invoked to understand the observed differences in the proper- 
ties of quasars, active galaxies, and BL Lac objects as due to 
the orientation of the relativistic beam and an obscuring torus 
with respect to the line of sight. Thus comparison of appar- 
ent jet speeds with optical classification is an effective test of 
these unified models and also serves to refine their parameters. 

4.2. Jet Kinematics 

The motion statistics collected for our sample address many 
physical questions related to relativistic jets. The process of 
their formation, that is, the initial acceleration and collima- 



tion of the flow, may be constrained by studying the speeds as 
a function of distance from the beginning of the visible jet (of- 
ten referred to as the "core"), and also by comparing the times 
when flares occur with the back-extrapolated epochs at which 
moving features appear to originate. Long-term multi-epoch 
observations can show whether there are accelerations or de- 
celerations farther down the jets, and whether radio features 
follow straight or curved trajectories. 

In order to discriminate between various jet models, and 
then to refine the relevant ones, it is important to establish 
whether jets exist as predefined channels, along which mul- 
tiple moving features can be seen, or whether instead suc- 
cessive components follow different trajectories, at either the 
same or at different speeds. Viable models of jet formation 
will also need to be able to reproduce the observed Lorentz 
factor distribution. The moving features may in fact be pat- 
terns, caused by the propagation of shocks rather than the 
flow of matter. This can be studied by comparing the Lorentz 
factors inferred from the motions to the Lorentz factors esti- 
mated by other means, such as variability, brightness temper- 
ature, and (possibly) the presence of X-ray and gamma-ray 
emission. By probing for correlations between the apparent 
velocities and other quantities such as the radio or X-ray lu- 
minosity, more can be learned about the physical parameters 
relevant for jet formation. Finally, studying the distribution 
of velocities as a function of optical host type is relevant to 
constraining unification models. 

Most of the jets in our study are well collimated and are 
unresolved transverse to their flow, although there often is 
significant curvature. In some radio galaxies with two sided 
structure, the appearance of the source is very frequency de- 
pendent, suggesting free-free absorption in a disk or torus, 
probably associated with the accretion disk surrounding the 
central engine. 

We found several jet features to have apparent negative ve- 
locities; that is, they appear to be approaching rather than re- 
ceding from the core. Most of these apparently negative ve- 
locities are consistent, within the errors, with no significant 
motion. Observations extending over a longer time frame are 
needed to determine if these inward motions are real. Appar- 
ent inward motion may be produced if there is a newly emerg- 
ing jet feature that is ejected from the core, and the combina- 
tion is not resolved by our beam. This would cause a shift 
in the measured position of the centroid and a corresponding 
decrease with time in the apparent separation of the core and 
other jet features. It is also possible that the true core is not 
seen, possibly due to absorption, and that the only observed 
features are parts of a jet which are moving with different ve- 
locities. If the furthest component is moving with a slower 
velocity than the one closest to the obscured core, then they 
will appear to be approaching each other. The apparent de- 
crease in component separation from the core could also be 
due to component motion away from the core along a highly 
curved jet which bends back toward the line of sight, so that 
the projected separation from the core appears to decrease 
with time. Finally, the moving features could be only pat- 
terns in the flow, some of which might even be moving back- 
ward. Istomin & Pariey ( 1996) have discussed an electron- 
positron jet model where an observer located close to the di- 
rection of motion will see backward-moving knots. We note 
that any model involving patterns must account for the the 
very large numbers of observed outward motions, as opposed 
to inward ones. For this reason, many of the simplest "moving 
marquee" models for superluminal motion have already been 
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Fig. 4. — Distribution of tlie apparent linear velocity in all sources with 
a quality code 'E' or 'G' and which have measured redshifts. The left hand 
column displays the distribution for all individual features which we have 
observed. Distributions in the center and right columns show only one feature 
per source, the brightest or the fastest respectively. Sources are divided by 
optical class in the second, third and fourth rows of the figure. 



ruled out. In the case of 0735+178, the most distant feature 
appears to be moving inward, but this is probably an artifact 
of the complex brightness distribution whose centroid shifts 
when the intensity distribution changes. 

4.2.1 . Velocity Distributions 

Figure 13 shows the distribution of the observed values of 
/3app = 1-58 X 1Q~^ fj, Da (I +z), where Da is the angular size 
distance to the radio source in megaparsecs and the angular 
velocity, /i, is in mas per year Figure|5]shows the same distri- 
butions for th ose j ets found in the representative sub-sample 
described in § 12.21 We include in these figures only those jet 
features which we have been able to measure with a quality 
code of 'E' or 'G' and which have measured redshifts. 

The central and right-hand columns of Figures|3and|5lshow 
the motions of only one feature per source, the brightest and 
the fastest respectively. The brightest feature in each jet is de- 
fined as the one with the largest peak flux density at the epoch 
for which the source image was presented in Papers I or II, 
although in nearly all cases this does not change over the pe- 
riod covered by these observations. These features generally 
have well-determined motions; however, in a few cases, the 
brightest jet feature does not have an 'E' or 'G' quality code. 



Fig. 5. — Same distributions as in Figure |3 except only those sources 
contained in our representative flux density-limited sub-sample are included. 



usually because it has a large diffuse structure. In these cases, 
we have substituted the brightest component which did have 
an 'E' or 'G' quality code. For the right-hand panel of Fig- 
ure |3 the fastest feature within a source is simply defined as 
the one with the largest (absolute) linear speed with a quality 
code of either 'E' or 'G'. 

K-S tests indicate that at a 95 percent confidence level, 
there is no difference in the distribution of the jet speeds, inde- 
pendent of whether we consider the brightest, fastest, or even 
all the jet features within a source. We consider the brightest 
feature of each source to be the most representative for our 
analysis, since in the case of sources with weak secondary 
features, observations having a limited sensitivity or dynamic 
range might only detect the brightest jet feature. 

The distributions shown in Figures |4] and |5] are also sub- 
divided according to optical class. We find that the velocity 
distributions for the quasars, radio galaxies, and BL Lacs are 
mostly concentrated in the same range (0 < /?app < 15). How- 
ever, the quasars have a tail ranging up to /3app ^ 34, while 
jets associated with active galaxies have a narrow range of 
velocities with /3app < 6. A K-S test confirms at the 98 per- 
cent confidence level that the quasars have a different speed 
distribution than the galaxies and BL Lac objects, while the 
distribution of speeds for the galaxies and BL Lacs are sta- 
tistically indistinguishable. Examination of Figure|4]and Fig- 
ure |5l suggests that the distribution of quasar velocities may 
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be bi-modal with a minimum near jSapp ^ 5. We have com- 
pared our observed distributions with one and two Gaussian 
distributions and find in each case a significantly lower value 
of reduced for the two component distribution for all of 
the quasar distributions. However, we are reluctant to quan- 
tify this further, as the observed distributions are clearly more 
complex than can be represented by two Gaussians. 

We have also compared the dispersion in the speed of dif- 
ferent features within each jet, with the dispersion of the av- 
erage jet speeds. We find that the dispersion of the average jet 
speeds among the 50 sources with two or more features with 
quality factors E or G is 6.3c, whereas in all but four jets (92 
percent) the dispersion in the speed of individual jet features 
is smaller than this value. This suggests that there is an under- 
lying flow, characteristic of each jet, that determines the speed 
of individual features with the jet. 

4.2.2. Apparent Speed and Apparent Luminosity 

Figure |6l shows the relation between /3app on apparent 
VLB A core luminosity at 15 GHz. There is a distinct up- 
per envelope to the distribution, which is very similar to that 
found a t 5 GHz in the Caltech Jodrell Flat-spectrum survey by 
IVerme ulen ( 1995). In particular, the low-luminosity sources 
all have slow apparent speeds. If we divide the sources at the 
median luminosity of Lig < 10^^ "^ WHz"', a K-S test yields 
less than a 1% probability that speed distributions for the 
brightest components of the high and low luminosity sources 
have the same parent distribution. We note that if the apparent 
jet speeds that we are measuring are simply random patterns, 
then no envelope w ould be expected in Figure El 

As discussed by iLister & MarscheJ J1997 ). the presence 
of this envelope does not necessarily imply that int rinsica lly 
faint jets have low intrinsic speeds. As discussed in § 14.3.11 in 
practice there is likely to be a spread in Lorentz factors. The 
lowest-luminosity sources in a flux-limited sample will tend 
to lie at low redshifts, where the co-moving volume element 
is small. Ignoring evolutionary effects and assuming a con- 
stant co-moving space density and a steep luminosity func- 
tion, these sources will therefore be representative of the most 
common jets in the parent population. The probability of any 
of them having both a high Lorentz factor and small viewing 
angle will be very small, especially if the parent population 
is dominated by low-Lorentz factor jets (i.e., n{j) oc 7"'-^ as 
found by Lister &Marscher 1997). This may therefore ex- 
plain why none of the low-luminosity sources in Figure|6lare 
highly superluminal. 

4.2.3. Non-Radial and Bent Trajectories 

For most of the sources in our sample we have determined 
only the radial speed and position angle of the motion with re- 
spect to the assumed core component. For those sources with 
sufficiently high-quality data, we also computed the full two- 
dimensional vector proper motion on the sky. This was ac- 
complished by fitting separately for the proper motion in both 
right ascension and declination. These results were combined 
to form a vector velocity, determined by the speed, /i, and the 
direction, 0, on the sky relative to the radio core. We are par- 
ticularly interested in comparing the direction of motion, (p, 
to the mean structural position angle of the jet, <i?>. 

Table|3]compiles the vector proper motions for all jet com- 
ponents in our sample which have at least five epochs of ob- 
servation and for which the vector velocity is of at least 5(t 
significance. These criteria guarantee that only the highest 
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Fig. 6. — Apparent speed of the fastest feature found in 110 jets as a 
function of the radio luminosity. 



quality vector motions are used in our comparison of velocity 
direction to structural position angle. Of the 60 component 
motions that meet these criteria, we find that 20 are signifi- 
cantly "non-radial", meaning the velocity direction, 0, differs 
by at least 3cr from the mean structural position angle, <d>. 
These non-radial motions are high-lighted in bold in Tabled 
Non-radial motions do not follow a direction which extrap- 
olates back to the jet origin and are by definition considered 
non-ballistic. The occurrence of such non-ballistic motions 
in approximately a third of the highest quality motions we ex- 
amined is striking. Assuming the observed velocities trace the 
underlying jet flow, this result indicates that bends in jet di- 
rection and/or jet collimation are common. The actual bends 
may be small, only a few degrees, but the observed bends can 
be large, of the order of 90°, since they appear amplified by 
projection. 

In FigureQ we show tracks in the RA - Dec plane of several 
sources that have clearly defined non-radial or bent trajecto- 
ries. In general, the trajectories are bent toward the next struc- 
ture down the jet, whether seen directly in our maps, or in pub- 
lished larger scale images. This suggests that jet features may 
trace out a flow in a pre-existing curved channel. 1226+023 
(3C273), 1219+285, 1532+016, and 2200+420 (BL Lac) are 
excellent examples of this trend. In 1548+056, the compo- 
nent motion appears to be transverse to the main jet direction 
which is known from lower resolution images, and is hinted 
at in our images, to lie to the north. However, the motion we 
observe is directed toward the small extension toward the east, 
and looking back toward the core, the jet ridge-line appears to 
curve into the component from the opposite direction. 

In 3C 390.3 (1845+797) the two most prominent features 
appear to be moving with a similar speed of about 2.5c, but 
along slightly different trajectories of 27° and 30°. Neither 
trajectory is aligned with the direction of the core or with the 
narrow 60" jet that points toward one of the distant hot spots 
in position angle 35° ( Alef et aL.1996») . This is a clear excep- 
tion to the trend identified above. 
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Fig. 7. — Selected images of sources with jets that show non-ballistic component motion. Measured component positions for each epoch are superimposed 
on the images along with the vector motion (solid lines) in the RA - Dec plane. Dashed lines represent the mean structural position angles, <i?>, for each 
component. 
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In general we do not have sufficient data to robustly de- 
tect changes in velocity, although the high incidence of non- 
ballistic motion implies that such changes must occur. How- 
ever, in BL Lac, there is evidence that the moving feature 
traces the curved ridge-line of the jet. Also, we have previ- 
ously reported an abrupt chang e in both speed and direction 
in 3C 279 that occurred in 1998 jHoman etal.l2003h . 

4.2.4. Two-sided Sources 

Seven sources in Table |2] display two-sided structure in 
which the jet flow appears bi-directional away from a central 
core. We have previously reported our results on the two- 
sided source NGC 1052 (Vermeulenet al. 2003a), where the 
component motion is close to the plane of the sky and mul- 
tiple components move in opposite directions from the core 
with a velocity, (3 ^ 0.26, which is only mildly relativistic. 
Plots of component position vs time for all of the two sided 
jets are shown in Figure|2] For most of these sources we were 
able to identify the central component at each epoch and so 
determine the motion of individual features away from the 
center. In a few cases, however, such as 1404H-286 (epochs 
1998.83 and 1999.55) the central core was not detected at 
one or more epochs, either because it was weak, or because 
the dynamic range was limited for that observation. In these 
cases we have interpolated from the positions of the outer fea- 
tures, to find the "virtual" center, and we have used that to 
determine the motion of the individual outer features. For the 
source 1323H-321 even that procedure was not possible since 
the structure appears to be two-sided, but we do not see the 
core at any epoch. Moreover, our imaging of this source is 
not completely satisfactory due to its large angular size and 
the limited interferometer spacings sampled by our data. 

4.2.5. Peaked Spectrum Sources 

Gigahertz Peaked Spectrum (GPS) sources are character- 
ized by their sharp low frequency cutoff and general absence 
of prominent large scale structure. There is no consensus 
as to whether the observed spectral cutoff in these sources 
is due to synchrotron self absorption or to free-free absorp- 
tion (l^haffer. Kellermann. & Cornwell 1999). Previous VLBl 
observations have suggested that GPS sources have a simple 
double structure morpho logy with little or no significant mo- 
tions (e.g., 0'De3[l99^. Thirteen of our sources are classi- 
fied as GPS sources of which we, so far, have observations 
of eight at three or more epochs. Of these, only 1345H-125 
shows superluminal motion. The mean velocity observed for 
the GPS sources in our sample is only (0.5 ± 1.5)c. None of 
the GPS sources have /3app > 1.5. 

The existence of a sharp peak in the spectra of these sources 
implies that the individual components each peak up at about 
the same frequency. This means that the different components 
probably have comparable Doppler shifts, and if the spectral 
cutoff is due to synchrotron self absorption they should have 
comparable brightness temperatures. Some GPS sources, 
such as 0552H-398 and 0642H-449 contain a single unresolved 
strong component, and so an interpretation in terms of syn- 
chrotron self absorption is required for any reasonable mag- 
netic field strength. However, for other GPS sources we note 
the very different appearance of the core and jet components; 
nevertheless, they have a common cutoff frequency. More- 
over, many of the peaked spectrum sources included in our 
study have sharply bent jet structure suggesting very different 
amounts of Doppler boosting and frequency shift of the self 



absorption peak. Free-free absorption probably plays an im- 
portant role in determining the spectra of these sources. Since 
many of the core-jet GPS sources show little or no motion, 
we suggest that possibly, the jet flow in these sources is non- 
relativistic, or that these sources are seen at very large viewing 
angles, rather than that the measured velocity is the advance 
speed of a young dou ble lobe radio galaxy as suggested by 
lOwsianik & Conwavl ( il998i) . Either of these scenarios would 
imply a small Doppler shift. GPS sources do not have any 
extended double lobe structure which may be a consequence 
of the slow jets which do not carry sufficient energy to form 
extended radio lobes. 

We note that in some cases, such as CTA 102 which had 
a peaked spectrum in the past, there are la rge variations in 
the total flux density jKovalev et al. 2002) which probably 
reflects significant relativistic boosting. Lister (2003) and 
Kovalev (2004) have argued that these sources are not bona 
fide GPS sources but "masquerading blazars" which often 
contain bright transient jet features that dominate their radio 
spectrum. 

4.2.6. Sources With Extended Double Structure 

As a result of our selection criteria, most of the sources 
included in our study are dominated by their flat spectrum 
compact structure; they are mostly identified with quasars. 
Some quasars, however, also have extended structure with 
relatively steep radio spectra, in addition to their flat spec- 
trum compact core. We have included a few of these lobe- 
dominated sources in our study, although they did not meet 
our selection criteria. Extended structure is more common 
among the active galaxies. The following sources in Table 121 
have significant double-lobe steep-spectrum extended struc- 
ture: NGC 315 (0055H-300), NGC 1052 (0238-084), 3C 111 
(0415H-379), 3C 120 (0430H-052), M87 (1228H-126), 3C 390.3 
(1845H-797), and Cygnus A (1945H-405). AH ai-e identified 
with galaxies containing a relatively strong AGN. All but 
NGC 1052 and 3C 120 are dominated by the extended steep 
spectrum double lobe structure even at 15 GHz. In the frame- 
work of standard relativistic beaming models, radio galaxies 
and lobe-dominated sources are presumed to lie close to the 
plane of the sky, and thus should show values of /J^pp ^ 1. 
However, 3C111, 3C 120, and 3C 390.3 each show rather 
large apparent motions, and it is not clear how they fit into 
standard unification models. 

These objects are all broad-line radio galaxies (BLRG), al- 
though 3C 120 has also been classified as a narrow-line ra- 
dio galaxy (NLRG). The NLRG are thought to lie at rather 
large angles to the line of sight, and polarimetry shows that 
in some cases they contain a central quasar and a broad 
line region that are hidden by a dusty torus. The NLRG 
do not show strong superluminal motion. 3C111, 3C 120, 
and 3C 390.3 do show superluminal motion and hence must 
be at small angles to the line of sight (less than about 11° 
for 3C 111). The BLRG are not well-understood, and they 
may have a wide range of intrinsic luminosity and orienta- 
tion( Dennett-Thorpe, Barthel, & van Bemmel 2000). 

On the other hand, the jets of other lobe-dominated sources 
appear to be sub-relativistic. The powerful radio galaxy 
Cygnus A contains twin jets pointed toward the distant hot 
spots. The approaching and receding jets appear to propagate 
with velocities of 0.7c and 0.2c respectively, and we find rel- 
atively slow speeds of (0.16±0.03)c in the parsec scale jet 
of NGC 315. In the case of the radio galaxy NGC 1052, we 
have reported a two-sided flow with sub-relativistic jet speeds 
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of only 0.26c JVermeulen et alJ^2003a^ on both sides of the 
core. At least in these radio galaxies, the intrinsic jet flow 
close to the core appears to have only moderate speeds with 
/3app <C 1. Thus, it appears that the observed distribution of 
speeds and/or degree of jet asymmetry cannot be interpreted 
simply in terms of the orientation of a twin relativistic jet with 
(3^1, unless the observed speeds do not reflect that actual jet 
flow speed, or the radiation from a high speed inner spine is 
not observed due to a more slowly moving outer sheath that 
dominates the emission at large angles from the jet direction. 
As we show in § l4.4.1l highlv relativistic jets are relatively rare 
among the general population of radio jets. 

4.2.7. Comparison With Other Velocity Studies 

Other recent VLBA monitoring observations (e.g. , 
iJorstad et al. »2001a'tf.'Bri tzen et alJ200UlHoman et alJ200U 
[Vermeulen et al. 2003b) made at other wavelengths with 
different resolution and with different sampling intervals 
are complementary to our observations and may be used 
to extend the range of size and time scales over which jet 
kinematics may be studied. 

Jorstad et al. (2001a b) have used the VLBA to study the 
motions in a sample of strong sources at 7 mm and 13 mm. 
These observations probe the source structure and motions on 
a scale about three times smaller than our 2 cm observations, 
but due to the decreased surface brightness at 7 mm, they are 
generally not able to trace the motions beyond a few milliarc- 
seconds from the core. From the comparison of the observed 
motions for the sources in common to our two studies, it is 
possible to trace the motions of individual components over a 
wider range of scales than is possible from either set of obser- 
vations alone. 

We have compared our velocities determined at 15 GHz 
with those of Jorstad et al. (2001a b). There appears to be 
little agreement in the individual s ource spe eds found by the 
two studies. On average, the Jors tad et alJ (12001 a b) speeds 
appear systematically higher, possibly because their observa- 
tions were carried out at shorter wavelengths and probed jet 
regions closer to the core. 

Brit zen et alJ d 19991 1200 ih and IVermeulen et'ai] ( l2003bh 
have discussed jet motions in a large sample of sources at 
6 cm, taken from the Caltech-Jodrell CJ surveys covering de- 
clinations greater than +35°. These observations have lower 
angular resolution than ours, and are thus more sensitive to 
the lower surfac e brightness structure located downstream. 
IVermeulen et alJ (l2003b,) quote a mean velocity for quasar 
and galaxy jets of 2.9c and 0.9c respectively. This appears 
be less than the values of (7.3±0.8)c and (1.7±0.8)c which 
we measure at 2 cm for the brig htest features in each source. 
On the other hand, Jorstad et al. ( 200 la), working primarily at 
0.7 and 1.3 cm, find systematically faster velocities in those 
sources where our samples overlap. These results suggest that 
there is a systematic decrease in /3app with increasing wave- 
length, probably because the observations at different wave- 
lengths, sample different parts of jet structure. In a separate 
paper (R. C. Vermeulen et al., in preparation) we will discuss, 
in more detail, the motions of those sources in common to the 
6 and 2 cm samples. 

4.3. Implications for Relativistic Beaming Models 

4.3.1. Ballistic Models 

The simplest model to consider is a pure ballistic model 
with a common flow velocity for all sources, with (3^ = ^ \ ; 



and a flow which is intrinsically symmetric. In this simple 
case, observations of the apparent velocities, and the ratio of 
flux densities of approaching and receding components, can, 
in principle, be used to solve uniquely for 7 and 9 and thus 
provide a test of the hypothesis that /3p = /3b. However even 
early VLBl data indicated that this simple mod el is no t tenable 
(Vermeulen & Cohen 1994; Lister & Marsch eill997h . 

In a flux density-limited sample, the combined effect of 
available solid angle and Doppler boostin g leads to a distribu- 
tion o f /9app with many values close to 7 ( Vermeule n & Cohei] 
119941) . On the other hand, if Doppler boosting is not an impor- 
tant selection mechanism, then in most cases /?app ^ 1 . Fig- 
ure |3l shows the expected distribution of apparent velocities 
for the case where 7p = 7b (see § 14. 1> along with the corre- 
sponding distribution for a randomly oriented sample. None 
of the distributions shown in Figure|4]or Figure|5]are consis- 
tent with either of these simple ballistic models as there is nei- 
ther evidence for the peak at /S^pp ^ 7 characteristic of the sim- 
ple Doppler boosted models, nor for the sharp low speed cut- 
off at /3app ^ 1 characteristic of models which ignore Doppler 
boosting. 

How do we reconcile the difference between the observed 
and predicted distributions? Either there must be spread in in- 
trinsic velocity (e.g.. Lister & Marscher 1997) or a difference 
between the bulk flow velocity and the pattern velocity, so that 
there is less Doppler bias in favor of observing beams which 
are or iented close to the critical angle ( Vermeulen & Cohe J 
In Figure|3j we plot a model with a spread of intrinsic 
velocity such that n(7) cx 7"' -^. This model has equal pattern 
and bulk velocities, and reproduces the general characteristics 
of the observed distribution of superluminal speeds. A more 
detailed analysis will be discussed by M. L. Lister et al. (in 
preparation). 

4.3.2. Randomly Oriented Samples 

Ba sed on e arlier observations of apparent velocity distribu- 
tions |Se^& Liang ( 1990) suggested that the simple model 
with no Doppler bias, provided an adequate fit to the ob- 
served distributions of /3app with only a slight adjustment 
needed which was satisfied by introducing an obscuring torus. 
However, comparison of our data with models which do not 
include Doppler boosting shows poor agreement, since, as 
shown in Figure |3l in the absence of Doppler boosting, most 
sources are expected to lie close to the plane of the sky where 
App ^ 1 ■ This appears to be inconsistent with the tail of the 
velocity distribution which we find extending toward larger 
apparent velocities. Following the discussion of Cohen ( 199(|) 
a detailed analysis of the observed velocities shows that the 
probability of the sources having been picked at random from 
a parent population that is isotropically distributed is less than 
10"^. We therefore conclude that the sources are preferen- 
tially aligned along the line of sight, as would be expected if 
Doppler boosting is in fact important. 

4.4. Comparison with Other Relativistic Velocity Indicators 

While the actual observation of superluminal motion 
in radio source jets remains the most direct way of es- 
tablishing the existence of relativistic motion, compari- 
son of the observed values of /3app with other veloc- 
ity indicators also provides important tests of beam- 
ing models and specifically the relation between pattern 
and bulk flow speeds. This include s flux density vari- 
abilitv iHughes et al. .1992: Lahteenmakr& Valtaoial 119991 
iLahteenmaki. Valtaoia. & Wiiklil999t) . maximum brightness 
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Fig. 9. — Histogram of Lorentz factors for 49 sources, calculated from 
/3app and Sy^ based on T'"' = 2 X lO'" K. 



Fig. 8. — Apparent velocity, ftpp, plotted against variability Doppler 
factor, <5var, fo r tlie fastest com ponent found in 49 sources calculated using 
the method of ILahteenmaki & V altaoia ( 1999) with an intrinsic biightness 
temperature TjJ"' = 2 X 10'" K. The solid line shows the expected locus of 
points for Lorentz factor values of 25. The dotted line represents the I/7 
cone where <5var = /3app. 



temperature (^Kellermann & Paulinv-Totfl 119691: iReadheadI 
[l994; Guiiosa & Dalv 1996), the ratio of core-to-extended 
radio luminosity (Orr& Browne 1982), the gamma-ray lu- 
minosity ( von Montiany et al. 1995; Hartmanetal. 1999; 
iMattoxet al. 1997, 2001), or observations of inverse Comp- 
ton X-rays (,Ghisellini et al. 1993). Relativistic boosting 
also affects the radio source counts and lumin osity functions 
(|Padovani & Urry 1992; Wall & .Iacksori ll99"% which provide 
a consistency check on beaming models. 

4.4. 1 . Relation Between Observed Velocity and Flux 
Density Variations 

Flux density changes are commonly seen in superluminal 
sources, and their short time scale is generally taken to imply 
high brightness temperatures. The variability time scale, and 
the time scale for apparent transverse motion, are both com- 
pressed due to the forward motion, and hence we might expect 
to see a relation between /5app and the flux density variabil- 
ity provided that /?b is related to /3p. ILahteenmaki & Valtaoial 
([1999) used variability data at 1 .3 cm and 8 mm from the Met- 
sahovi Observatory to calculate a variability Doppler factor, 
(5var, assuming an intrinsic brightness temperature characteris- 
tic of a self-absorbed synchrotron source in which the particle 
and magnetic energy are in equilibrium. We have recalculated 
their values using the cosmology given in §nand for different 
values of intrinsic brightness temperature. Since (5var varies in- 
versely as the cube root of the assumed intrinsic temperature 
in the synchrotron plasma, T^^^, (5vai is relatively insensitive to 
the assumed value of T^'"'. 

In Figure |8] we plot /3app against (5var for the 49 sources in 
common to the Metsahovi and VLBA samples. We calcu- 
late (5var = [7j^v3Y7"jjnt]'/3 assuming a spectral index of zero and 
where T^^"'^ is the apparent brightness temperature calculated 
from the variability time scale, by assuming that it is limited 
by the size of the source divided by the speed of light. Five 
sources which only have components located at a bend in the 
jet were excluded, as they probably reflect a standing shock 
wave, or perhaps a stationary location in a heUcal jet where 



several values of 7j^"" 



the flow is closest to the line of sight, and hence is boosted 
most strongly. In either case the measured velocity is a poor 
indicator of the flow velocity, and not relevant in discussing 
relativistic effects. The sample we use contains 5 active galax- 
ies, 14 BL Lacs, and 30 quasars. Most points lie inside the 
"1 /7 cone" (/3aDD = i^var) as they should for a flux desnity lim- 
ited sample ( Vermeul en & Cohenll994t) . 

Cohen et al.. C2003») have compared Figure |8] with the re- 
sults of a simulation generated by randomly picking a flux- 
limited sample from an isotropically distributed population 
with power law distributions of luminosity and 7 and for 
over the range 4 x 10'' to 1 x 10'' K. 
2 X 10'" K gives the best fit between the measured 
and simulated data. However, although there does appear 
to be an upper limit to /3app that is close to the expected 
locus for components with 7 = 25, the detailed distribution 
is not well-matched to that expected from the simulation 
dLister & Marscheill997h . 

Calculations of the variability Doppler factor using values 
of T™^ closer to the inverse Compton limit, ^ 10'^ K, lead 
to distributions on the /3app-<^var plane that are very different 
from those expected from the simulations. We conclude that 
T™^ is perhaps an order of magnitude below the inverse Comp- 
ton limit. 

For this application, we have used the fastest feature for 
each source, on the grounds that these velocities should be 
more representative of the true flow velocities. Slower- 
moving components, especially those located at a bend in the 
jet, may be dominated by standing shock waves. Forward 
shock waves might also exist, and trying to understand their 
role is a goal of our study. Other geometries have been sug- 
gested for the jet, including a fast "spine" which we would 
preferentially see, surrounded by a slower shell. In this case, 
the spine would probably also control the flux variations, so 
that using the fastest (spine) velocity for /3app is appropriate. 

We have also examined values of (5var calculated from the 
UMRAO data base at our wavelength of 2 cm, and noted a 
large dispersion between Doppler factors deduced from these 
data and the shorter wavelength Metsahovi data. Thus, the ro- 
bustness of the Doppler factors calculated in this way appears 
to be uncertain. The Michigan data cover a longer time span, 
but are at longer wavelengths where individual outbursts ap- 
pear to overlap in time. 

The Lorentz factor, 7, in a superluminal jet is important, be- 
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cause it is intrinsic to the jet, whereas /3app and (5var are the ob- 
servables that depend on the jet orientation. The distribution 
of 7 may give information on the physics of the collimation 
region. From the values of /3app and 5 shown in Figure Owe 
calculate the actual Lorentz factors, shown as a histogram in 
Figure |9l We believe that the apparent deficit in the first bin 
of Figure |9] for the quasars, may be a selection effect or just 
due to small number statistics. To be fully consistent with the 
analysis above, and the choice 7jJ'" = 2 x 10^" K, the Lorentz 
factors should have a power-law distribution. However, the 
number of objects is too small, and the errors are too uncer- 
tain, to make a meaningful comparison. A small Lorentz fac- 
tor implies a low velocity and small flux density variations, 
and the latter especially is less likely to be measurable. How- 
ever, the calculated Lorentz factors are useful in showing that 
there must be a wide range of 7 in the superluminal sources. 
The galaxies all have rather low Lorentz factors, while the 
quasars have a broad distribution up to 7 sa 30. 

4.4.2. Gamma-Ray Sources 

Many of the sources included in our study have been 
cataloged as strong gamma-ray sources according to mea- 
surements made by the EGRET detector on board the 
Comp t on Gamma-Ra y Obs e rvatory ( von Montianv et al. 
[1995; Hartman et al.' 1999; Mattox et al. i 19971. 1200 It 
ISow ards-Emmerd. Romani, & Michelson 2003). It is gener- 
ally thought that the gamma-ray emission occurs deep within 
the relativistic jet. Arguments based on size limits deduced 
from time variability and the cross section for pair production 
suggest that the gamma-ray emission, like the radio emission, 
is Doppler boosted ( Dermer & Schlickeiser 1994). In fact, 
the gamma-rays may be even more strongly beamed than 
the radio emission since the former generally have a steeper 
spectral index, a, and the flux density boostin g varies as 
j2-Q continuous jets. Also, Permed il995h has shown 
that if the bulk of the gamma-rays are produced by external 
Compton scattering off photons associated with the accretion 
disk, the resulting gamma-ray emission will be boosted by an 
additional factor of J'"". 

If gamma-ray loud AGN do indeed have systematically 
high Doppler factors, then we might expect to see a differ- 
ent apparent speed distribution for them than for AGN that 
have not been detected in gamma-rays. The situation is com- 
plicated by the possibility that the gamma-ray loud jets may 
be seen inside the critical angle for maximizing superlumi- 
nal motion, (I/7 radians), and therefore might have slow 
apparent projected speeds. However, Monte Carlo simula- 
tions based on a simple linear relationship between rad io and 
gamma-ray luminosity (e.g.,|Lister & Marscher 1999: lListed 
[1999) confirm that in a flux density limited radio sample, 
AGN detected by EGRET should have typically higher speeds 
than t hose that w ere not detected in gamma-rays. 

Jorstad et alJ J2001ah recently measured the apparent 
speeds of 33 EGRET-detected AGN. They found a mean value 
of 16c for the fastest component in each source and con- 
cluded that the gamma-ray sources have larger Lorentz fac- 
tors than th e general population of radio sources. However, 
LTorstad eFaL (.2001a) did not have a non-gamma-ray control 
sample observed in the same way and at the same wavelength 
with which to compare their results. 

We have classified our sample into EGRET 
a nd non- EGRET sources based on the list of 
attoxetal. (2001.) and recently mod ified by 
Sowards-Emmerd. Romani. & MichelsonI J2003I) . These 
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Fig. 10. — Histogram of the brightest component speed in EGRET and 
Non-£Gfi£r detected sources for our representative flux density limited sam- 
ple. 



authors c lassify gamma - ray so urces from the third EGRET 
catalog jHartman et alJ Il999h as "highly probable" and 
"plausible" AGN identifications, based on Bayesian statistics 
and their proximity to bright flat-spectrum radio sources. 
In Figure ^| we show the distribution of measured speeds 
for the EGRET detected and non-detected sources. For the 
brightest jet component found in the 18 EGRET sources 
in our representative sample we found a median speed of 
(8.0 ± 1.6)c compared with a value of (3.9 ± l.l)c for the 
53 sources with no EGRET detections. A K-S test suggests 
that the difference appears significant at the 90 percent level. 
For the purpose of this analysis we have included the two 
"plausible" EGRET sources, 2230H-114 and 1156H-295 as 
detections. Re-classifying them as non-detections sources 
had no effect on our results, and neither did excluding the 
sources with negative velocities from our analysis. 

These results are consistent with the idea that the radio 
emission from gamma-ray quasars is indeed more strongly 
beamed than for the whole radio quasar population. However, 
our samples are incomplete and may therefore be biased. 

4.5. The Angular Velocity Redshift Relation 

FigureFTTIshows the measured values of angular velocity, 
vs redshift, for the fastest 'E' or 'G' rated component found 
in the sources in Table|2l The line represents /imax, the fastest 
proper motion a source can display if it has 7 = 25. The vari- 
ables are the observables, uncontaminated by modeling, and 
hence are of value in showing directly that (a) at all redshifts, 
the observed velocities are not clustered near the maximum 
value as expected from the simple ballistic models with a sin- 
gle Lorentz factor for all sources, (b) low values of ijl are seen 
at all z, and (c) high values of fi are seen only at low z- This is 
true for galaxies, BE Lacs, and quasars and appears inconsis- 
tent with non-cosmological interpretation of quasar redshifts 
as proposed by Burbidge ( 2004). 

An early version of this plot iCohen et alJI 19881) was used 
to show that the standard paradigm for superluminal motion 
provided a crude upper bound to the points in Figure ^2 and 
that the maximum value of /3app in 32 sources was about 13 
(//o = 70, ^0 = 0.5). The /imax line in Figure [TTI is. similarly, 
a crude upper bound to the measured points, and shows that 
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Fig. 1 1. — fj,-z diagram showing tlie distribution of angular velocity for 
the fastest component of the 110 sources which has a quality factor of 'G' or 
'E'. The solid line is the maximum value of /3app for 7 = 25. 



Papp.max ~ 25. Notc that the cosmology used in this paper 
makes /3app somewhat larger than for a cosmology with Jim = 1 
and = 0. 

The data in iCohen etaUfml were compiled at assorted 
frequencies, mostly below 15 GHz, and corresponded to mea- 
surements made at different distances from the core which 
ma y expla in the difference in Papp.mux- As we have discussed 
in § l4.2.7l there is some evidence that observations at high fre- 
quencies give higher values of fj, th an lower frequencies. 

IVermeulen & Cohe n' ('1994) and Lister & Marscher (1997^ 
have shown that, if Doppler boosting is important, then 
even rather small samples of superluminal sources will show 
/3app,max near the maximum value of the distribution of 7, 7n,ax- 
Hence we expect that 7max ^ 25 for the sample shown in Fig- 
ure ^2 This is marginally consistent with Figure |9j reflect- 
i ng that Lorentz factors calculated according to the method in 
§ l4.4.1l have rather large uncertainties. 

5. SUMMARY 

We have studied the kinematics of a large well-defined sam- 
ple of 1 10 quasar and active galaxy jets and find a distribution 
of apparent velocities typically between and 14c but ranging 
up to about 34c. There is evidence for a characteristic velocity 
in each jet which may represent the true plasma flow velocity. 
We have found that quasar jets generally have larger apparent 
velocities than jets associated with BL Lac objects and active 
galaxies, although the distributions overlap. Our measured 
values of /S^pp are consistent with the Doppler factors, (5var, 
calculated from time variability and a parent population hav- 
ing a steep power law distribution of intrinsic Lorentz factors 
extending down to moderate velocities and an intrinsic bright- 
ness temperature close to 2 x 10'" K. This is close to the value 
expected if the particle and magnetic energy densities in the 
jet^are comparable (Readhead 1994; Singal & Gopal-Krishna 
Il985l) 

In approximately one-third of the well studied jets, we find 
evidence for non-ballistic motion; that is the flow is not along 
the direction away from the core. In most of the jets, we find 
no deviation from a constant speed; but in a few sources, we 



do see evidence for changes in speed and direction of indi- 
vidual features. Mostly, the flow appears to lie along the di- 
rection of the local jet orientation. However, in some cases 
the flow has a significant non-radial component, which points 
toward more distant parts of the jet. This suggests that there 
is a continuous flow along a pre-existing channel. Contrary to 
the assumption of the simple unified models, in some jets, the 
intrinsic flow appears to be with speeds much less than c. 

Observations made at higher frequencies sample jet fea- 
tures located closer to the core, and they typically show larger 
apparent velocities than we observe at 15 GHz, while lower 
frequency observations show yet smaller speeds. Sources 
with stable GPS spectra show little or no motion; the jet flow 
in these sources may be non-relativistic or lie in the plane of 
the sky with a correspondingly small Doppler shift. 

We find that jets of quasars which have been observed as 
strong gamma-ray sources have marginally higher speeds than 
those which are not gamma-ray sources. This is consistent 
with models where the gamma-ray sources have more highly 
relativistic jets and are aligned closer to the line of sight. 
However, our analysis is limited by both small number statis- 
tics, and the uncertainties in the ever changing analysis of the 
EGRET catalogs. Also, with the limited range of flux density 
observed by EGRET, there is no well defined class of gamma- 
ray loud and gamma-ray quiet sources analogous to the radio 
loud and radio quiet classifications. More sensitive observa- 
tions with the next generation of gamma-ray observatories, 
such as GLAST, combined with jet speed data for our com- 
plete radio sample of 133 radio sources, should be very useful 
for investigating gamma-ray production mechanisms in AGN 
jets and relating the gamma-ray properties to the observed jet 
outflow. 

E. Ros et al. (in preparation) have extended these obser- 
vations and analysis through 2001 and 2002. M. L. Lister 
et al. (in preparation) have defined a more complete sam- 
ple of 133 sources and are continuing the observations of 
these sources including linear and circular polarization (D. C. 
Homan et al., in preparation). Observations with this new 
sample will allow a more robust comparison with models, 
a better estimate of the distribution of intrinsic Lorentz fac- 
tors, and a start to understanding the evolution of jet magnetic 
fields. 
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TABLE 1 
General Source Properties 



lAU 




Sub-sample 


Optical 




Radio 


Radio 




EGRET 


Name 


Alias 


Member 


Class 


Redshift 


Spectrum 


Morphology 




ID 


(I) 


(2) 




(4) 


("51 


(6) 


(7) 




^91 


0003-066 


NRAO 005 


Y 


B 


0.347 


Flat 


CJ 


2.84 


N 


0007+106 


in Zw 2, Mrk 1501 


Y 


G 


0.0893 


Flat 


N 


1.29 


N 


0016+731 




Y 


Q 


1.781 


Flat 


CJ 


0.98 


N 


0026+346 




N 


G" 


0.517 


Flat 


CJ 


0.65 


N 


0035+413 




N 


Q 


1.353 


Flat 


a 


0.53 


N 


0048-097 




Y 


B 




Flat 


N 


1.65 


N 


0055+300 


NGC315 


N 


G 


0.0165 


Flat 


CJ 


0.82 


N 


0106+013 




Y 


Q 


2.107 


Flat 


CJ 


2.97 


N 


0112-017 




N 


Q 


1.365 


Flat 


a 


0.82 


N 


0119+041 




N 


Q 


0.637 


Flat 


CJ 


1.28 


N 


0133+476 


DA 55 


Y 


Q 


0.859 


Flat 


CJ 


4.95 


N 


0149+218 




N 


Q 


1.320 


Flat 


CJ 


1.42 


N 


0153+744 




N 


Q 


2.341 


Flat 


CJ 


0.37 


N 


0202+149 


4C +15.05 


Y 


Q' 


0.405 


Flat 


CJ 


2.29 


Y 


0202+319 




Y 


Q 


1.466 


Flat 


CJ 


2.27 


N 


0212+735 




Y 


Q 


2.367 


Flat 


CJ 


2.69 


N 


0218+357'' 




N 


B" 


0.96 


Flat 


CJ 


0.70 


N 


0234+285 


CTD 20 


Y 


Q 


1.213 


Flat 


CJ 


4.04 


Y 


0235+164 




Y 


B*" 


0.940 


Flat 


N 


1.62 


Y 


0238-084 


NGC 1052 


Y 


G 


0.00490 


Flat 


2S 


2.48 


N 


0316+162' 


CTA21 


N 


G 




Peaked 


CJ 


0.34 


N 


0316+413= 


3C 84 


Y 


G 


0.01756 


Flat 


CJ 


10.59 


N 


0333+321 


NRAO 140 


Y 


Q 


1.263 


Flat 


CJ 


2.24 


N 


0336-019 


CTA26 


Y 


Q 


0.852 


Flat 


CJ 


3.44 


Y 


0355+508 


NRAO 150 


N 


u 




Flat 


CJ 


7.11 


N 


0415+379 


3C 111 


Y 


G 


0.0485 


Steep 


CJ,LD 


5.98 


N 


0420-014 




Y 


Q 


0.915 


Flat 


CJ 


10.60 


Y 


0430+052 


3C 120 


Y 


G 


0.033 


Flat 


CJ 


4.41 


N 


0440-003 


NRAO 190 


N 


Q 


0.844 


Flat 


a 


1.20 


Y 


0454+844 




N 


B 




Flat 


a 


0.34 


N 


0458-020 




Y 


Q 


2.286 


Flat 


CJ 


2.33 


Y 


0528+134 




Y 


Q 


2.070 


Flat 


CJ 


7.95 


Y 


0552+398 


DA 193 


Y 


Q 


2.363 


Peaked 


N 


5.02 


N 


0602+673 




N 


Q 


1.970 


Flat 


N 


0.97 


N 


0605-085 




Y 


Q 


0.872 


Flat 


CJ 


2.80 


N 


0607-157 




Y 


Q 


0.324 


Flat 


CJ 


7.26 


N 


0615+820 




N 


Q 


0.710 


Flat 


N 


0.48 


N 


0642+449 


OH 471 


Y 


Q 


3.408 


Peaked 


CJ 


4.31 


N 


0707+476 




N 


Q 


1.292 


Flat 


CJ 


0.63 


N 


0710+439 




N 


G 


0.518 


Peaked 


2S 


0.61 


N 


0716+714 




Y 


B 




Flat 


CJ 


1.25 


Y 


0727-115 




Y 




1.591 


Flat 


CJ 


5.30 


N 


0735+178 




Y 


B 




Flat 


CJ 


1.64 


Y 


0736+017 




Y 


Q 


0.191 


Flat 


CJ 


2.58 


N 


0738+313 


Ol 363 


Y 


Q 


0.630 


Flat 


CJ 


2.87 


N 


0742+103 




Y 


G'' 




Peaked 


a 


1.50 


N 


0745+241 




N 




0.409 


Flat 


CJ 


0.95 


N 


0748+126 




Y 


Q 


0.889 


Flat 


CJ 


3.25 


N 


0754+100 




Y 


B 


0.266 


Flat 


CJ 


1.83 


N 


0804+499 




Y 


Q 


1.432 


Flat 


CJ 


1.14 


N 


0808+019 




Y 


B 




Flat 


CJ 


1.34 


N 


0814+425 




Y 


B 


0.245 


Flat 


CJ 


1.28 


N 


0823+033 




Y 


B 


0.506 


Flat 


CJ 


1.40 


N 


0829+046 




Y 


B 


0.180 


Flat 


CJ 


1.35 


Y 


0850+581 




N 


Q 


1.322 


Flat 


CJ 


0.61 


N 


0851+202 


OJ287 


Y 


B 


0.306 


Flat 


CJ 


4.12 


Y 


0859-140 




N 


Q 


1.339 


Steep 


CJ, CSS 


1.58 


N 


0906+015 


4C +01.24 


Y 


Q 


1.018 


Flat 


CJ 


2.73 


N 


0917+449 




N 


Q 


2.180 


Flat 


CJ 


1.43 


P 


0923+392 


4C +39.25 


Y 


Q 


0.698 


Flat 


CJ 


12.68 


N 


0945+408 




Y 


Q 


1.252 


Flat 


CJ 


1.59 


N 


0953+254 




N 


Q 


0.712 


Flat 


CJ 


1.31 


N 


1012+232 


4C +23.24 


N 


Q 


0.565 


Flat 


a 


1.16 


N 


1015+359 




N 


Q 


1.226 


Flat 


CJ 


0.82 


N 


1049+215 


4C +21.28 


N 


Q 


1.300 


Flat 


a 


1.45 


N 


1055+018 


4C +01.28 


Y 


Q 


0.888 


Flat 


CJ 


5.30 


N 


1055+201 


4C +20.24 


N 


Q 


1.110 


Flat 


CJ 


0.38 


N 


1101+384 


Mrk 421 


N 


B 


0.031 


Flat 


CJ 


0.53 


Y 


1127-145 




Y 


Q 


1.187 


Flat 


CJ 


3.39 


N 


1128+385 




N 


Q 


1.733 


Flat 


CJ 


1.13 


N 


1155+251 




N 


G» 


0.202 


Flat 


CJ 


0.24 


N 


1156+295 


4C +29.45 


Y 


Q 


0.729 


Flat 


CJ 


2.36 


P 


1219+285 


W Comae 


N 


B 


0.102 


Flat 


CJ 


0.60 


P 
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lAU 




Sub-sample 


Optical 




Radio 


Radio 


SvLBl 


EGRI 


Name 


Alias 


Member 


Class 


Keasniit 


Spectrum 


Morphology 


(Jy) 


ID 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


1226+023 


3C 273 


Y 


Q 


0.158 


Flat 


CJ 


41.40 


Y 


1228+126 


M87 


Y 


G 


0.0044 


Steep 


CJ, LD 


2.97 


N 


1253-055 


3C 279 


Y 


Q 


0.538 


Flat 


CJ 


24.89 


Y 


1302-102 




N 


Q 


0.286 


Flat 


CJ 


0.71 


N 


1308+326 




Y 


Q 


0.997 


Flat 


CJ 


3.98 


N 


1323+321 


4C +32.44 


N 


G" 


0.370 


Peaked 


2S 


0.65 


N 


1 328+254' 


3C 287 


N 


n 

V 


1.055 


c ^ 


n CSS 


0.09 


N 


1 'iTQj-'imf 


-iC OCA 
jL- zoO 


M 
iN 




A 9/1 A 
U.o40 


Steep 






iM 


i — iZ / 




V 
I 


V 




Fill 
rial 




0.0 / 


V 
1 


1 '^A'^J.1 0'^ 
i J'+J+lZJ 


+1^ +iZ.jU 


IN 




A 1 T 1 
U. izi 


Peaked 


7C 
ZO 


A QQ 


IN 


1 ACiA^I'kf, 


00 TAG AAG 
KJ\l zUo, iVilK OOo 


IN 




A A'77 


Peaked 


T^ 
ZO 


1 on 

l.zU 


IN 


1 Al 'l-l-l 
14-1 J+l J J 




V 
1 


13 


A lAl 
U.Z4- / 


Plat 


Zo 


1 10 


iM 






IN 




1 AQl 


TJlnt 

r^iai 


IN 


A A1 
U.Oi 


r 


1 ^riQ n^^ 




IN 




1 1 Ql 


Steep 


CJ, Coo 


A 71 


w 

i>l 


1 ^ 1 n ASQ 




V 

I 




A lAA 


rial 


CJ 


O'X 

z.yj 


"V 
I 


1 ^ T j_n 1 A 

i J Jz+UlO 




IN 




1 /ITA 


rial 


CJ 


A 7A 

u. /o 


IN 


i J40+UZ / 




v 

1 


/-V 


A /I 1 T 


nai 


CJ 


7 Ql. 
Z.OJ 


IN 


1 s j_n^A 

1 j'+o+UjO 


/ir^ j_A^ A/1 

+UJ.O+ 


V 

I 




1 /m 

1 .+ZZ 


plot 

rial 


CJ 


7 Q7 

z.yz 


W 
i>l 


1 AAA j_ 1 AA 
1 OUO+ i uo 




V 




1 TTA 
1 .ZZO 


plot 

rial 


CJ 




"V 




DA 406 


V 


n 

V 


1 401 


Flat 




JZ 


Y 


1633+382 


4C +38.41 


Y 





1.807 


Flat 


CJ 


4.28 


Y 




NR AO 1 2 


Y 


V 


1.666 


Flat 


N 


1.61 


N 


1 64 1 +3QQ 


3C 345 


Y 


n 


0.594 


Flat 


CJ 


8.73 


N 


1 642+690 


4C +69.21 


N 


V 


0.751 


Flat 


CJ 


1.34 


N 


1652+398 


Mrk 501 


N 


B 


0.033 


Flat 


CJ 


0.90 


N 


1655+077 




Y 


Q 


0.621 


Flat 


CJ 


2.09 


N 


1656+053 




N 


Q 


0.879 


Flat 


CJ 


0.69 


N 


1656+477 




N 


Q 


1.622 


Flat 


CJ 


1.14 


N 


1730-130 


NRAO 530 


Y 


Q 


0.902 


Flat 


CJ 


10.97 


Y 


1749+096 


4C +09.57 


Y 


Bb 


0.320 


Flat 


CJ 


5.58 


N 


1749+701 




N 


B 


0.770 


Flat 


CJ 


0.79 


N 


1758+388 




Y 


Q 


2.092 


Peaked 


CJ 


1.75 


N 


1800+440 




Y 


Q 


0.663 


Flat 


CJ 


1.50 


N 


1 803+784 




Y 


B'' 


0.680 


Flat 


CJ 


2.54 


N 


1 807+698 


3C 371 


N 




0.050 


Flat 


CJ 


1.38 


N 


ioZJ+JDo 


/H^ j-^A TT 
+JO.Z / 


v 
1 


D 


A AAl 


plot 

rial 


CJ 


Z.J I 


IN 


ioZo+H-o / 




V 




u.oyz 


Steep 


V.J, \^oo 


Z.UI 


iM 




lOA 1 


IN 




A A^7 


Steep 


CJ, 


A /17 


iM 




3C 395 


IN 








V.00 


L.J J 


N 


1921—293 




N 


o 


0.352 


Plat' 


CJ 


14.39 


N 




AC +73 1 8 


Y 


n 


0.303 


Flat 


CJ 


3.92 


N 


1957+405 




Y 


G 


0.056 




2S LD 


1.68 


N 






Y 




V 


1 736 


Flat 


CJ 


2 79 


N 


2007+777 




N 




0.342 


Flat 


CJ 


1.16 


N 


2021+317 


AC +31 56 


Y 


u 




Flat 


CJ 


2.16 


N 


9091+614 

ZVJZ 1 TU 1 't 




Y 


V 


227 


Flat 


CJ 


2 73 


N 


21 13+293 




N 


n 


1.514 


Flat 


CJ 


0.94 


N 


2128—123 




Y 


n 


0.501 


Flat 


CJ 


3.18 


N 


0111 AT 1 
Z i J i— Uzi 


^1.. — UZ.O i 


v 
1 


Rb 


i .ZOJ 


plot 

rial 


CJ 


01 

Z.Zl 


iM 


1 "5/1 j_AA/1 
Z i j't+UU'^ 




V 

I 




1 .yjz 


Peaked 


CJ 


A 'lA 


IN 


z i jO+i+i 


WA. 101 


V 

I 




7 /177 
Z.4Z / 


plot 

rial 


CJ 


7'^ 
Z. / J 


iM 


2144+Uyz 




JN 




1 111 


plot 

rial 


CJ 


A CI 


IN 


2145+067 


4C +06.69 


Y 





0.999 


Flat 


CJ 


10.37 


N 


2200+420 


BL Lac 


Y 


B 


0.069 


Flat 


CJ 


5.67 


Y 


2201+315 


4C +31.63 


Y 


Q 


0.298 


Flat 


CJ 


3.28 


N 


2209+236 




Y 


Q" 


1.125 


Flat 


CJ 


1.62 


Y 


2223-052 


3C446 


Y 


Q 


1.404 


Flat 


CJ 


6.57 


N 


2230+114 


CTA 102 


Y 


Q 


1.037 


Flat 


CJ 


4.86 


P 


2234+282 


CTD 135 


N 


Q 


0.795 


Flat 


CJ 


1.44 


N 


2243-123 




Y 


Q 


0.630 


Flat 


CJ 


2.56 


N 


2251+158 


3C 454.3 


Y 


Q 


0.859 


Flat 


CJ 


12.08 


Y 


2345-167 




Y 


Q 


0.576 


Flat 


CJ 


2.54 


N 



^ Source not in the lVeron-Cettv & Verorll2001l) catalog. 

^ Source classified as quasar in the''Veron-Cettv & Verorf i"2001|) catalog. 

Source classified as galaxy in the Veron-Cettv & Veron 1 20011) catalog. 
^ The source 0218+357 is a gravitationally lensed AGN which we do not include in our statistical analysis. ^^^^^^^^^^^ 

^ The complex source 0316+413 (3C 84) satisfies our selection criteria, but was observed as part of a separate monitoring program IWalker et aU20QCl) . It is included in our statistical 
analysis for completeness. 

^ Source has large or diffuse structures that are not adequately sampled by our observations. 
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TABLE 1 — Continued 



lAU 




Sub-sample 


Optical 




Radio 


Radio 


SvLBI 


EGRET 


Name 


Alias 


Member 


Class 


Redshift 


Spectrum 


Morphology 


(Jy) 


ID 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(V) 


(8) 


(9) 



Note. — C olumns are as follows: (1) lAU Name (B1950.0); (2) Other Name; (3) Indicator for the flux density limited sub-sample (see §1^; (4) Optical classification according to 
the lV^ron-Cett v & Veron i 2001 ) catalog, where Q - quasar, B - BL Lac object, G = active galaxy, a nd U - unid entified; (5) Redshift from Veron-Cettv & Veron ■ 2001 ) {the red shifts 
for 00 26+346, 0727-1 15 and 1 155+251 are from'Zensus et al."2002', that of 0202+149 is from Stickel eTalTTggl 0218+357 is from Lawrence 1996, 0754+100 is from Carangelo etall 
|2063 and that of 2209+236 is from Sowards-Emmerd. Romani, & Michelson 2003); (6) Description of radio spectrum (see §|2|; (7) Radio morphology classification, where CJ= 
core-jet, 2S = two-sided jet, N = naked core, LD = lobe -dominated, and CSS = comp act-steep spectium object; (8) Strongest 2 cm total cleaned VLBA flux density at any epoch, in 
Jy; (9) EGRET gamma-ray source identification according to lMattox et all 1200 ll and lSowards-Emmerd. Romani. & MichelsoJ 120031) . where Y = highly probable identification, P = 
probable identification, N = no identification. 

TABLE 2 

Proper motions. 



Ohiprt 


Comp. 


N 


< R> 


< 1? > 




ft 

Mapp 




Tlnztinr't rTimnrinpnt 


T? n ti n CT 








(4) ^ 




( mas/vrl 


(ycsrs) 






(1) 


(2) 


(3) 




(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


0003-066* 


Bt 


5 


1.1 


-72 


-0.05 ±0.09 


-1.1 ±2.0 




Y 


G 




C 


4 


3.0 


-77 


0.01 ±0.02 


0.3 ±0.4 




Y 


E 




D 


5 


6.2 


-78 


-0.02 ±0.06 


-0.3 ±1.3 




Y 


E 


0016+731* 


Bt 


4 


0.8 


-50 


0.074 ±0.003 


5.7 ±0.2 


1987.2 ±0.5 


N 


G 


0026+346 


B 


5 


10.4 


-124 


0.03 ±0.23 


0.9 ±7.0 




N 


P 




ct 


5 


12.2 


-129 


0.15±0.12 


4.7 ±3.8 




Y 


G 




E 


5 


14.7 


56 


-0.16±0.12 


-4.8±3.6 




Y 


G 




F 


4 


15.8 


56 


-0.02 ±0.02 


-0.6 ±0.6 




Y 


E 




G 


5 


18.2 


59 


0.02 ±0.05 


0.5 ± 1.5 




Y 


E 


0035+413 


Bt 


5 


1.5 


104 


0.10±0.02 


6.3 ± 1.5 


1982.6 ±3.5 


Y 


E 


0055+300 


Bt 


5 


0.9 


-54 


0.14±0.03 


0.16±0.03 


1991.0±1.1 


N 


G 




C 


6 


3.7 


-51 


0.00 ±0.04 


0.00 ± 0.05 




Y 


E 




D 


6 


6.0 


-51 


0.03 ±0.03 


0.03 ±0.03 




Y 


E 


0106+013* 


Bt 


4 


2.4 


-119 


0.28 ±0.04 


23.3 ±3.6 


1988.3 ±1.4 


Y 


E 




C 


5 


2.1 


-113 


0.27 ±0.03 


22.3 ±2.7 


1990.9 ±1.0 


N 


G 


0112-017 


Bt 


5 


1.2 


118 


0.02 ±0.03 


1.0± 1.8 




Y 


E 


0119+041 


Bt 


5 


0.4 


124 


0.01 ±0.04 


0.5 ± 1.6 




N 


G 


0133+476* 


Bt 


5 


2.7 


143 


0.04 ±0.01 


2.0 ±0.7 




Y 


E 


0149+218 


B 


3 


2.8 


-13 


0.03±0.18 


1.9±11.7 




Y 


P 




C 


6 


5.0 


-22 


0.29 ±0.03 


18.4± 1.9 


1980.8 ±1.8 


Y 


E 




Dt 


6 


8.8 


-12 


0.27 ±0.07 


16.8±4.6 


1964.9 ±9.0 


Y 


E 


0153+744 


Bt 


6 


10.4 


152 


0.01 ±0.02 


1.0± 1.5 




Y 


E 


0202+149* 


B 


3 


0.8 


-67 


0.25 ±0.06 


6.2±1.6 


1994.7 ±1.0 


N 


F 




Dt 


4 


5.1 


-59 


-0.0I±0.0I 


-0.2 ±0.3 




Y 


E 


0202+319* 


Bt 


3 


7.2 


-2 


0.06 ±0.03 


3.9 ±2.0 




Y 


G 


0212+735* 


Bt 


4 


2.4 


104 


0.08 ±0.01 


7.0± 1.0 


1966.1 ±4.5 


Y 


E 




C 


4 


14.0 


92 


-0.01 ±0.04 


-1.1 ±3.8 




Y 


E 


(19 34.+? 8 S* 


Bt 


5 


3 5 


— 13 




13 5 ± 2 9 


1QS7 0-1-3 4 


Y 






E08 


7 






Afi 4- ID 


n 9^ 4- 03 




Y 


E 




E17 








7S -1- 19 


n 9^ 4- 04 


1991 4di 9 


Y 


E 




CjZXJ ' 


/ 






U. /V It u. iz 


U.zO It U.U4 




V 

1 


1 ■ 
1 . 




E23 








7/^4-0 19 


9^-1-0 04 


1 QQ4 1 -1- n S 


Y 


E 




W07 


7 






0.82±0.13 


0.27 ±0.04 


1984.7 ±2.0 


Y 


E 




WIO 


7 






0.72±0.11 


0.24 ±0.04 


1991.2±1.0 


Y 


E 




W14 


7 






0.75±0.12 


0.25 ±0.04 


1993.4 ±0.7 


Y 


E 




W16 


6 






0.80±0.12 


0.27 ±0.04 


1995.4 ±0.4 


Y 


E 


0316+413*<> 


Bt 


3 






0.19±0.05 


0.2±0.1 


1960.0 ±5.0 


Y 


E 


0333+321* 


Bt 


9 


2.7 


123 


0.18±0.01 


11.0±0.5 


1982.7 ±0.7 


Y 


E 




C 


9 


4.9 


126 


0.20 ±0.03 


12.0± 1.7 


1972.5 ±3.5 


Y 


E 




D 


3 


7.2 


128 


0.40 ±0.07 


24.5 ±4.5 


1980.0 ±3.3 


Y 


G 




E 


4 


11.1 


131 


-0.08 ±0.13 


-4.7 ±8.1 




N 


F 


0336-019* 


Bt 


5 


1.4 


67 


0.22 ±0.04 


10.1 ±2.0 


1991. 8± 1.4 


Y 


E 




C 


4 


1.6 


46 


0.21 ±0.08 


9.7 ±3.9 




Y 


G 


0355+508 


Bt 


5 


3.2 


46 


0.09 ±0.03 






Y 


E 


0415+379* 


B 


4 


2.3 


65 


1.29±0.32 


4.2± 1.0 


1994.5 ±0.5 


Y 


P 




C 


3 


4.2 


70 


0.90 ±0.53 


2.9± 1.7 




Y 


P 




Dt 


6 


3.3 


66 


1.52 ±0.05 


4.9 ±0.2 


1996.7 ±0.1 


Y 


E 


0420-014* 


Bt 


5 


1.6 


-163 


0.03 ±0.07 


1.5±3.6 




N 


G 




C 


5 


3.1 


170 


0.29±0.14 


14.2 ±6.6 




Y 


G 


0430+052* 


Bt 


10 


4.7 


-111 


1.77 ±0.06 


3.9±0.1 


1995.0±0.1 


Y 


E 




C 


8 


6.1 


-106 


1.80±0.19 


4.0 ±0.4 


1993.7±0.4 


Y 


E 




D 


3 


6.1 


-110 


1.36±0.35 


3.0±0.8 


1991. 8± 1.2 


Y 


P 




G 


9 


3.0 


-112 


1.59±0.11 


3.5 ±0.2 


1996.5 ±0.2 


Y 


E 




I 


3 


1.7 


-125 


1.51 ±0.44 


3.4±1.0 


1999.0 ±0.4 


Y 


P 




H 


5 


2.6 


-118 


2.08 ±0.24 


4.6 ±0.5 


1998.2 ±0.2 


Y 


E 


0440-003 


Bt 


4 


1.1 


-120 


0.03 ±0.03 


1.1±1.4 




Y 


E 
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Object 


Comp. 


N 


<ff > 


<i?> 






^0 


Distinct component 


Rating 








(mas) 


(deg) 


(mas/yr) 


(years) 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


0454+844 


Bt 


5 


1.3 


162 


-0.10 ±0.05 






Y 


E 


0458-020* 


Bt 


3 


1.5 


-43 


0.10 ±0.04 


8.9± 3.3 




Y 


G 




c 


3 


2.8 


-51 


0.11 ±0.05 


9.4±4.1 




Y 


G 




Bt 


3 


4 


23 


n 077 4- f) 009 


6.4 ± 0.2 


1 Q4.fl 0-\- \ 


Y 


G 


UDUJ — KJoJ 




3 


1 6 


123 


10^091 


4-. J HI y .o 




N 


p 




rt 


Q 


J.O 


1 1 n 


1 S ^ 09 


O.O ^ U.O 


1 Q77 Q -1- 9 


V 




OdO/— Ij / 




4 


2.4 


88 


0.1 / ±U.U4 


3.5 ± 0.8 


iyo4.1 ± 3.5 


Y 


E 


0642+449 


Tit 

B ' 


3 


0.2 


92 


0.01 ± 0.02 


0.9 ± 2.0 




N 


F 


07074476 


B ' 


4 


3.5 


5 


— 0.U4 it O.Ui 


—2.8 ±1.7 




Y 


E 




C 


3 


O.J 


24 


— U. 1 1 zh U.Uj 


—7.1 ±3.3 




Y 


KJ 


0710+439 


BT 


4 


13.3 


3 


0.01 ±0.13 


0.2 ± 3.9 




Y 


G 




C 


4 


14.4 


3 


-0.05 ± 0.07 


-1.4 ±2.0 




N 


G 




D 


4 


15.8 


4 


0.00 ±0.02 


0.0 ±0.5 




Y 


E 




G 


4 


24.4 


1 


0.00 ± 0.01 


0.0 ± 0.3 




Y 


E 


U710+714* 


nt 


4 


2.4 


14 


0.65 ± 0.04 




1994.2 ± O.J 


N 


G 




rj 

a 


5 


1 Q 
1 .O 


1 1 
1 1 


n ^ 1 _t- A A/1 
U.J 1 ±: U.U4 




1 QO^ ^ _U A 1 

lyyj.j it U.J 




r* 


0/2/— 1 15 


B ' 


3 


2.7 


-52 


0.44 ± 0.01 


31.2±0.6 


1990. / ± 0.1 


Y 


G 


0735+178* 


BT 


4 


0.7 


72 


0.17 ±0.02 




1992.5 ±0.7 


N 


G 




C 


3 


1.2 


59 


f\ ^ A 1 /^^ 1 ^ 

0.64 ±0.10 




1994.2 ±0.5 


N 


P 




L) 


c 

J 


4.U 


04 


A 1 Q _L A A1 

— U. lo it U.Uj 






I 


c, 




R 


O 




7 A 
/U 


A nq _L A A9 
u.vj Hi u.uz 


11 -1- A 9 
i i . J ^ U.Z 


1 QQ'^ S -1- A 1 

ly vj. J It u. i 


V 


c 




J) 




10.0 


—88 




4 8 ± 8 


1972 7 ±4 1 


Y 






Ft 


J 


Z. 1 


OO 


A 8 S -1- A 1 A 
U.O J ^ U. iU 


1 A ^ -1- 1 
iU.J It i . J 




V 


c 


0/3o+315 


15 ' 


8 


2.0 


179 


A ri/i _L A A1 

U.UO lb O.UI 


2.3 ± 0.4 


1 nci _L A 

lyjj.o ± O.O 


Y 


E 






7 


'X 1 


1 7A 
1 /u 


A A7 -1- A 1 A 


AA-'X Si. 
Z.*+ III J.O 




IN 


r 




Rt 


D 




—J 


A (M. _l_ A A1 
U.UJ It U.Ui 






IN 


n 
\j 




Q 




0.7 


—38 


_A 09 -k 09 






N 


G 


U / +Ji-Z+1 


rt 


o 
o 


T 1 


—Oj 


A 19 -U A A^ 
U. JZ It U.Uj 


7 Q _|_ 1 1 

/.y it i .J 


1 QCQ 7 _l_ 1 A 

lyoo.z It 1 .0 


I 


C 




[) 


g 


6 


-64 


A 001 A- 099 


0.1 ± 0.6 




N 


G 


U /4o+lzO 


Tit 


6 


1 .7 


120 


/^^ 1 o _L Am 

U. Iz it U.UI 


5.9 ± 0.5 


1 nOQ Z _L 1 1 
lyOJ.J ±1.1 


Y 


E 






Q 


4 7 


1 1 1 
111 


A 97A -1- A AAA 
U.Z /'4 ^ U,UU4- 


1 'I 1 -1- A 9 
i J, i It U.Z 


1 QSA 8 -k A 9 

lyou.o dz U.Z 


V 


KJ 


U /j4+lUU 


T3t 

15 ' 


5 


A A 
4.U 


15 


A IPt _L A 1 O 

U. /U it U. Iz 


1 1.9 ± 2.1 


1 nni 1 _L 1 A 
lyyi.j ±1.0 


Y 


E 




L- 





A Q 


1 / 


A A^ _L A <\1 
U.Uj it U.U/ 


A Q J- 1 9 
U.O it l.Z 




JN 


KJ 


0804+499 


B ' 


4 


0.8 


134 


0.13 ± 0.06 


8.5 ±4.2 




N 


G 


0808+019* 


B' 


3 


0.4 


-175 


0.11 ±0.01 




1993.1 ±0.5 


N 


F 




C 


3 


4.8 


-173 


-0.20 ±0.1 8 






Y 


P 


0814+425* 


BT 


6 


0.9 


90 


0.18±0.02 


2.9 ±0.3 


1993.2 ±0.5 


N 


G 




C 


5 


1.7 


129 


0.32 ±0.11 


4.9 ± 1.7 




N 


F 


1 (Til* 


C 


6 


2.4 


14 


f\ AO 1 /^^ ("1 /I 

0.48 ± 0.04 


14.4 ±1.1 


1993.3 ±0.4 


N 


G 






6 


1.1 


32 


0.31 ±0.06 


9.5 ± 1.8 


1994.8 ±0.8 


N 


G 




E 


3 


0.7 


32 


0.41 ±0.01 


12.3 ±0.4 


1998.2 ±0.1 


N 


F 


0829+046 


B 


4 


1 .3 


66 


0.30 ± 0.04 


3.5 ± 0.4 


1994.0 ±0.6 


N 


G 




C' 


5 


3.0 


66 


0.13 ±0.06 


1.5 ±0.7 




Y 


E 


0850+581 


B 


3 


0.8 


174 


0.05 ± 0. 1 1 


3.0±7.1 




N 


P 






4 


2.5 


152 


0.04 ±0.01 


2.5 ±0.5 


1935.6 ± 12.5 


Y 


E 




D 


4 


7.0 


151 


0.20 ± 0.07 


12.7 ±4.1 


1963.5 ± 11.3 


Y 


E 


0851+202 


C 


5 


1 .0 


—92 


0.52 ± 0.04 


9.9 ± 0.9 


1 on A 1 \ r\ 

1994. / ±0.2 


N 


G 




D' 


5 


0.9 


-102 


0.37 ±0.06 


7.1 ± 1.1 


1994.9 ±0.4 


N 


G 




r; 
E, 


4 


A A 


11/; 

— 110 


U.J 1 it U.Uz 


O.U it U.4 


1 QO^ "7 _L A 1 

lyyj. / it u.i 


IN 


KJ 




R 


-J 


1 1 


1 /^A 


U. lO It U.UM- 


1 A 1 -1- 9 S 
lU, 1 It Z.o 


1 QQA 'X -\- ^ Q 

lyyu. J It i .y 


IN 


P 

r 




rt 


4 


j.Z 


— 1 JJ 


A A1 -LA A 1 
U.Uj it U.UI 


1 .y it U.O 


1 QQQ Q _L ■29 9 

looy.o it j2.2 


I 


C, 




D 


3 


7.0 


158 


A O/i _L A 1 ^ 

0.zO±0.1j 


16.3 ± 9.4 




Y 


F 


OyOD+015 


rt 


4 


0.9 


38 


0.10±0.0/ 


5.1 ± 3.9 




N 


G 




D 


4 


1 A 

1 .4 


40 


U. Iz it U.U/ 


/: A \ '-> c 
0.4 it J.J 




N 


r~< 
KJ 




E 


4 


2.8 


45 


AT) J_ A (\n 

U.ZZ ± U.U/ 


1 1.8 ± 3.9 


lyoO.j ±4.3 


N 


/~i 

Cj 


0917+449 


B ' 


6 


1.2 


178 


0.07 ± 0.01 


5.8 ±0.9 


1979.8 ±2.6 


Y 


E 


092j+j92 


Tit 


6 


1.9 


—78 


0.0/ ± O.Oj 


2.9 ±1.0 




Y 


E 






7 


2.5 


—79 


A CI 1 _i_ A CIQ 

U.UI it U.Uj 


A T _L 1 O 
U.Z it 1 .Z 




N 


KJ 


0945+408 


BT 


5 


0.1 


-83 


0.03 ±0.04 


2.0 ±2.4 




N 


G 




C 


3 


8.3 


119 


0.37 ±0.20 


22.5 ± 12.1 




N 


P 


Uyjj+zD4 


13 




1 .Z 


1 AQ 

— lUo 


1 1 1 A A1 
U.J 1 it U.UI 


IT A J- (\ A 
LZA ± U.4 


1 om /I _L A 1 

iyy2.4 it u. 1 


IN 


r* 




rt 


J 


1 A 


1 9Q 


A AQ -1- A A1 


J. J It U.'^- 


lyoo. i It i . J 


NT 
IN 


KJ 




D 


3 


1.6 


-104 


0.301 ±0.001 


12.09 ±0.04 


1994.17±0.02 


N 


F 




E 


3 


0.4 


-124 


0.06 ±0.02 


2.4 ±0.7 


1992.0 ±2.2 


N 


F 


1012+232 


Bt 


5 


2.2 


110 


0.27 ±0.02 


9.0 ±0.6 


1989.9 ±0.5 


Y 


E 


1015+359 


Bt 


6 


1.3 


-172 


0.20 ±0.04 


12.4 ±2.2 


1991.0±1.2 


N 


G 




C 


5 


3.1 


-178 


0.04 ±0.04 


2.6 ±2.2 




Y 


E 




D 


7 


4.7 


-178 


0.16 ±0.04 


9.9 ±2.5 


1969.4±7.1 


Y 


E 


1049+215 


Bt 


6 


1.1 


106 


0.01 ±0.07 


0.9 ±4.4 




N 


G 




C 


6 


7.4 


112 


0.14 ±0.04 


8.8 ±2.6 


1944.7 ±16.1 


Y 


E 
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TABLE 2 — Continued 



Object 


Comp. 


iV 


<R> 


<i?> 


fJ'r 




to 


Distinct component 


Rating 








(mas) 


(deg) 


(mas/yr) 


(years) 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(V) 


(8) 


(9) 


(10) 


1055+018* 


BT 


5 


1.7 


-49 


0.04 ± 0.03 


1.7± 1.6 




Y 


E 


1055+201 


Bt 


7 


1.2 


-10 


0.02 ± 0.02 


0.9± 1.1 




N 


G 




C 


7 


3.0 


-5 


0.18 ±0.08 


10.0 ±4.3 




Y 


E 


1101+384 


B 


5 


2.2 


-40 


0.40 ±0.15 


0.8 ±0.3 




N 


F 




C 


7 


1.4 


-39 


0.22 ± 0.02 


0.46 ± 0.04 


1991.5 ±0.5 


N 


G 




Dt 


6 


0.8 


-30 


0.17 ±0.03 


0.4 ±0.1 


1993.8 ±0.7 


N 


G 




E 


3 


0.3 





0.07 ± 0.03 


0.1 ±0.1 




N 


F 




F 


3 


5.7 


-44 


-0.02 ± 0.30 


0.0 ± 0.6 




N 


P 


1127-145* 


B 


6 


2.3 


82 


0.05 ±0.06 


3.0±3.5 




Y 


E 




ct 


7 


4.1 


83 


0.09 ±0.01 


5.4±0.8 


1952.5 ±6.9 


Y 


E 


1128+385 


Bt 


5 


0.3 


-168 


0.01 ±0.01 


1.1±0.5 




N 


G 




C 


5 


1.0 


-155 


0.004 ± 0.008 


0.3 ±0.6 




N 


G 


1155+251 


B 


3 


10.2 


8 


0.02 ±0.17 


0.2 ±2.2 




N 


P 




C 


3 


10.2 


20 


-0.17 ±0.14 


-2.2 ± 1.9 




N 


P 




D 


3 


8.0 


14 


0.00 ± 


0.0 ±4.3 




N 


P 




ET 


3 


3.4 


21 


0.03 lb 0.03 


0.3 ±0.4 




Y 


G 




F 


3 


4.5 


-83 


-0.06 ± 0.20 


-0.7 ±2.5 




N 


P 




G 


3 


3.0 


-51 


0.01 ±0.02 


0.1 ±0.3 




N 


F 




I 


3 


1.2 


-24 


-0.03 ± 0.04 


-0.4 ±0.6 




N 


F 


1156+295* 


Bt 


8 


3.8 


-2 


0.22 ±0.08 


8.9 ±3.4 




Y 


G 


1219+285 


B 


6 


0.8 


101 


0.08 ±0.08 


0.5 ±0.6 




N 


F 




C 


7 


2.5 


107 


0.02 ±0.08 


0.1 ±0.5 




Y 


E 




Dt 


7 


11.0 


107 


0.48 ±0.03 


3.2±0.2 


1974.5 ±1.4 


Y 


E 


1226+023* 


Bt 


11 


6.1 


-113 


1.05 ±0.04 


10.8 ±0.5 


1992.1 ±0.3 


Y 


E 




C 


7 


9.6 


-120 


1.36 ±0.11 


14.1 ±1.1 


1991.3 ±0.6 


Y 


E 




D 


8 


10.8 


—119 


O.O J ± 0.20 


8.6 ±2.0 


1 no A 1 _i_ 1 1 


Y 


P 




E 


10 


3.3 


— llo 


1 T7 _L_ n 
1.2/ ±0.0/ 


13.1 ±0.0 


iyyj.o±o.2 


Y 


E 




f 


z: 
D 


10. / 


1 T2 

— IZj 


0. /y ± 0. ID 


0.2it l.D 


ly /4.0 it 4.0 


v 


ii 




Ij 


4 


2.3 


1 


0.41 ± o.jy 


4.2 ± 4.U 




Y 


P 




I 


4 


1 .4 


—122 


O.oo ± 0.32 


9.1 ± 3.3 




Y 


P 


1 OOQ_i_l OA* 


r 1 

D 


y 


o.u 


— 10 


0.14 ± 0.0/ 


n r\A _i_ n m 
U.U4 it U.U2 




rs 


r* 
\j 




r>t 
C 


10 


1.5 


—78 


— 0.01 diO.Oo 


n An _i_ A AO 
0.00 ± 0.02 




N 


G 


1253-055 


B 1 


10 


3.3 


-116 


0.28 ±0.01 


8.9 ±0.4 


1986.2 ±0.6 


Y 


E 


1302-102 


BT 


3 


1.5 


30 


0.31 ±0.05 


5.6 ±0.9 


1992.3 ±0.8 


N 


F 


1308+326* 


Bt 


3 


0.8 


-80 


0.313 ±0.002 


16.3±0.1 


1995.08 ±0.01 


Y 


G 


1323+321 


C 


3 


56.4 


-51 


0.04 ±0.03 


0.9 ±0.8 




N 


F 




D 


3 


58.0 


-50 


0.142 ±0.001 


3.25 ±0.02 


1588.5 ±1.4 


N 


F 




E 


3 


14.6 


-46 


-0.86 ±0.35 


-19.7 ±8.0 




N 


P 




F 


3 


10.8 


-41 


-0.54 ±0.28 


-12.3 ±6.5 




N 


P 




G 


3 


4.8 


-35 


-0.10 ±0.12 


-2.2 ±2.8 




N 


P 




it 


3 


55.0 


-50 


-0.15 ±0.03 


-3.4 ±0.7 




N 


F 


1334-127* 


Bt 


3 


2.5 


156 


0.05 ±0.03 


1.7 ±1.0 




N 


F 


1345+125 


B 


5 


1.1 


168 


0.09 ± 0.02 


0.7 ±0.2 


1986.8 ±3.2 


Y 


E 




C 


5 


2.0 


164 


0.15 ±0.02 


1.2±0.1 


1985.7 ±1.5 


Y 


E 




D 


5 


7.1 


171 


-0.07 ±0.09 


-0.6 ± 0.7 




Y 


G 




Et 


5 


9.8 


162 


0.03 ±0.01 


0.2 ±0.1 




Y 


E 




F 


5 


20.6 


156 


0.16 ±0.23 


1.3±1.9 




Y 


P 




G 


5 


44.5 


155 


-0.01 ±0.01 


-0.06 ±0.05 




Y 


E 




1 


5 


52.4 


161 


0.02 ± 0.05 


0.1 ±0.4 




Y 


E 




J 


5 


66.1 


177 


0.03 ±0.08 


0.2 ±0.6 




Y 


E 


1404+286 


BT 


10 


3.3 


53 


0.04 ± 0.02 


0.2 ±0.1 




Y 


E 




C 


9 


2.0 


62 


0.06 ± 0.03 


0.3 ±0.2 




Y 


E 




D 


10 


3.3 


—127 


O.UO ± 0.02 


0.0 ±0.1 




Y 


E 




J::. 


1 n 
lU 


/I A 
4.4 


—111 


U.Uj zt U.U2 


n 1 1 A 1 
U. 1 ± U. 1 




V 


ii 


1 /I 1 1 1 1 1 C* 

141 j+1 j5 


T5t 


6 


2.1 


—113 


0.25 ± 0.08 


3.9 ±1.2 


lyoy.o ± 1.1 


Y 


E 








j.a 


1 1 1 

— Ill 


n -17 _L n IS 

\}.J/ ZiZ\}. iD 


J.o ZLl 2.4 




V 


p 




F 


3 


6.6 


-109 


4i + 15 


7.1 ±2.4 


1984.9 ±4.9 


Y 


F 




Bt 


7 


1.7 


85 


V./ . 1 V./ ZD V./ . LJZ. 


6.2± 1.2 




Y 






Q 


5 


4 9 


83 


V./. J J ZD V./. JV./ 


31.4 ± 18.1 




Y 


p 


ijiu— uoy 




J 


Z.a 


—20 


U.oj it U.Uj 


1 Q O ^ 1 A 

lo.y ± 1 .u 


iyy4.D it 0.2 


V 


ii 






3 


2.3 


—32 


O.J / lb O.Oi 


12.0 ± (J. / 


iyyD.o±o.2 


N 


F 


1532+016 


B ' 


6 


1.1 


130 


0.21 ± 0.01 


13.7±0.9 


1992.6 ±0.4 


Y 


E 




u 


/: 
O 


0.0 


1 /IT 
14J 


U.Uo it U.U3 


A A \ A f\ 

4.1 ± i.y 




JN 


r* 
\j 


1546+027* 


Bt 


3 


1.2 


175 


0.05 ± 0.03 


1.3±0.8 




N 


F 


1548+056* 


Bt 


5 


4.4 


10 


0.052 ± 0.004 


3.5 ±0.3 


1912.6 ±7.0 


Y 


E 


1606+106* 


Bt 


5 


0.9 


-58 


0.30 ± 0.02 


18.4±1.4 


1994.8 ±0.3 


N 


G 




C 


4 


1.2 


-53 


0.38 ±0.03 


23.2 ±2.1 


1994.0±0.3 


N 


G 




D 


3 


1.5 


-52 


0.36 ±0.02 


21.6±1.3 


1991.9±0.2 


N 


F 




E 


3 


2.1 


-40 


0.50 ± 0.03 


30.1 ±1.6 


1991.9±0.2 


N 


F 




F 


3 


0.6 


-61 


0.196 ±0.001 


11.8±0.1 


1996. 16 ±0.00 


N 


F 


1611+343* 


Bt 


5 


2.9 


173 


0.06 ±0.04 


4.2 ±2.6 




Y 


E 
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TABLE 2 — Continued 



Object 


Comp. 


iV 


<R> 


<i?> 






to 


Distinct component 


Ratinj 








(mas) 


(deg) 


(mas/yr) 


(years) 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


1633+382* 


ot 
B' 


8 


1.5 


-86 


0.15 ±0.03 


1 1.5 ± 2.0 


1987.5 ± 1.8 


Y 


E 






7 

1 


z.o 


87 
— oZ 




7 ^ 7 ft 
/ . J Z.O 




L\ 


ri 




Rt 


Q 


1 7 




n /I Q _i_ f ^ 


1 7 ^ n Q 
1 / .u ± u.y 


LyyD.y zt U.Z 


V 
1 


c 

d 








n Q 


—101 


V^. J / ZL yiXiD 


1 2 8 ± 1 


1 QQ'^ A -1- 7 

177 J ."4- ^ U.Z 


Y 






Rt 

ID ' 


J 


7 7 


1 A7 
— iOZ 


{\ Q g _i_ n H/l 
vJ. JcS zt 


1 ft n 4- 1 s 


1 QQl 0-1-0 7 

lyyi.y it u. / 


v 
1 


c 
d 




J) 


4 


9 9 


—166 




—1 6 ± 2 3 




Y 




1 652+398 


B 


5 


2.7 


142 


0.06 ±0.26 


0.1 ±0.6 




Y 


p 








7.5 


109 


0.09 ± 0.05 


0.2 ±0.1 




Y 


E 


1655+077* 


B 


5 


3.7 


-38 


0.43 ±0.03 


15. 3± 1.0 


1988.9 ±0.6 


Y 


E 




ct 


5 


8 


—44 




2 4 ± 9 




Y 




iOJO+Uja 


Rt 


+ 


n 7 
u. / 


77 
/Z 


n no -U n m 


/I /I 4- n s 

4.4 It U.o 


1 QQO O -U 1 A 

lyyu.u it 1.+ 


IN 


KJ 


1 /^S(^-ud77 


Rt 




1 1 
i . i 


n 


n -I- n 


7 8 4- n 
Z.O ni J.u 




Ay 


ri 




Q 




7 


— 19 


VJ.UO ZL V^.UO 


4 5 ± 4 6 




Y 




1 l'X(\ 1 "^A* 
1 / JU l^v 


Rt 




^ n 

J.U 


1 7 


i\ A9. +-019 


77 7 4- < ft 
Zj.J It J.O 


1 QS7 S -1- 7 
lyo / . J It Z.J 


V 

1 


ri 
yj 


1 749+096* 


B 


4 


4.6 


29 


0.06 ± 0.07 


1.2± 1.3 




Y 


E 




ct 


4 


7 


34 


n 1 S ± 0'^ 

\} . ZC "J . ytJ 


3 1 ± 0.5 


IQQ'l C _|_ A O 
lyyj.j ^ u.o 


N 


Q 


1 /4y+/ui 




6 


2.3 


—66 


U.Ui ± U.Uj 


0.4 ±1.1 




Y 


E 








Q 7 


<c 

—JO 


— U.Uo4 it U.UUJ 


7 ft 4- n 1 

—J.O it U. i 




JN 


1 ■ 
r 






a 


1 ^ 
i .J 


7ft 


n 08 -!- n Of^ 

U.Uo It U.UD 


7 zl 4- 7 ft 
J. 4 IE Z.O 




"M 
i> 


r 


1 / j(S+j(Si5 


Rt 


4 


n A 
U.4 


— yo 


U.UUz it U.Ui J 


n 7 4- 1 1 
U.Z it i . i 




JN 


Lj 


1800+440* 


BT 


5 


2.6 


-163 


0.49 ±0.05 


18.5 ±2.0 


1992.5 ±0.6 


Y 


E 




C 


4 


4.4 


-155 


0.56±0.14 


21.2±5.3 


1989.1 ±2.0 


Y 


G 


1803+784^ 


J. 

BT 


6 


1.4 


-94 


-0.01 ±0.01 


-0.5 ± 0.4 




Y 


E 


1807+698 


Bt 


7 


2.5 


-96 


0.01 ±0.04 


0.0 ±0.1 




Y 


E 




C 


8 


5.2 


-94 


0.12 ±0.04 


0.4 ±0.1 




Y 


E 




D 


3 


6.8 


-97 


0.85 ± 0.47 


2.9± 1.6 




Y 


P 




E 


4 


0.6 


—97 


i\ 7^ ^ f\ 
Vj.ZD ± U.Ufi 


0.9 ± 0.3 


1 nn£ n 4_ n o 
1996.9 ±0.6 


N 


F 


1823+568* 


BT 


7 


1.7 


-158 


-0.003 ± 0.040 


-0.1 ± 1.5 




Y 


E 




C 


7 


6.7 


-161 


0.09 ± 0.03 


3.4± 1.1 


1923.0 ±23.4 


Y 


E 


l82s+4s7 


B 


3 


1.2 


—20 


0.z3 ± 0.03 


9.2± 1.3 


1993.5 ± 0. / 


N 


F 




\^ 


a 

J 


J.J 


7^ 
— Zj 


n ^ R ^ n 1 ^ 

\).jo niy). LJ 


14.0 zn J.O 




V 


r 




nt 


J 


y.o 


77 
— JZ 


U.zD ± U.Ui 


in 1 n A 
iU. i ± U.4 


iVOi . J zt 1 .D 


V 

1 


Li 


lo4D+/y / 


B 


9 


5.5 


—41 


U. J4 ± U.Uj 


2.1 ±0.1 


1 OCT o _L n 
iVo /.z ± U.O 


Y 


E 




CT 


10 


1.1 


-36 


0.60 ±0.01 


2.3±0.1 


1985.7 ±0.3 


Y 


E 


1901+319 


fit 


7 


0.8 


122 


0.03 ± 0.02 


0.9 ±0.7 




N 


G 




c 


7 


15.8 


120 


0.00 ±0.01 


0.0 ± 0.5 




Y 


E 


1921-293 


B 


4 


5.7 


24 


0.19 ±0.06 


4.2± 1.3 


1968.9 ±8.9 


Y 


E 




C 


4 


2.3 


-1 


0.17 ±0.06 


3.7± 1.2 


1984.7 ±4.5 


N 


G 




Dt 


4 


0.5 


-29 


0.01 ±0.01 


0.2 ±0.2 




N 


G 


1928+738* 


Bt 


8 


2.1 


165 


0.29 ±0.03 


5.6 ±0.6 


1990.3 ±0.7 


Y 


E 




C 


6 


3.6 


161 


0.30 ±0.08 


5.8 ± 1.5 


1984.7±3.1 


Y 


E 




E 


5 


0.8 


156 


0.12±0.06 


2.3± 1.1 




Y 


E 


1957+405* 


B 


5 


6.0 


-78 


0.12±0.02 


0.5±0.1 


1950.6 ±8.0 


N 


G 




D 


8 


3.0 


-79 


0.18 ±0.02 


0.7 ±0.1 


1981.4± 1.6 


N 


G 




E 


8 


2.2 


-80 


0.19 ±0.03 


0.7 ±0.1 


1986.0± 1.6 


N 


G 




F 


6 


1.7 


-79 


0.22 ±0.03 


0.8±0.1 


1 r\f\r\ 1111 

1990.3 ± 1.1 


N 


G 




GT 


8 


1.2 


-79 


0.16 ±0.03 


0.6 ±0.1 


1990.0 ± 1.6 


N 


G 




I 


7 


0.8 


-79 


0.10 ±0.02 


0.4 ±0.1 


1990.6 ± 1.6 


N 


G 




J 


5 


0.5 


—123 


A AC _1_ A AO 

U.U5 ± 0.02 


0.2 ±0.1 


1 non a _L 7 1 
19s9.6±3.1 


N 


G 




Pv 


7 
/ 


D.D 


1 no 


A no -4- n 1 A 


n n 4- n *\ 




V 

I 


ri 




L 


g 


3.0 


115 


A n< _|_ A AO 


0.2 ±0.1 




Y 


E 


2001+403* 


Bt 


5 


3.2 


120 


16±004 


12.3 ±3.0 


1978 2±4 8 


Y 


E 


2007+776 


Bt 


g 


1.4 


-92 


00 + 03 


0.0 ±0.6 




Y 


E 




D 


6 


6.8 


-91 


-0.04 ±0.23 


-0.8 ±4.9 




Y 


P 


2021+317* 


Bt 


5 


3.4 


166 


0.14 ±0.05 






Y 


E 




C 


3 


1.0 


-172 


0.02 ±0.02 






N 


F 


2021+614* 


B 


6 


3.9 


34 


0.04 ±0.04 


0.6 ±0.5 




Y 


E 




Ct 


6 


3.1 


-149 


0.05 ±0.04 


0.8 ±0.6 




Y 


E 


2113+293 


Bt 


5 


0.3 


177 


0.02 ±0.01 


1.4 ±0.3 


1985.1 ±3.2 


N 


G 


2128-123* 


Bt 


3 


3.7 


-143 


-0.07 ±0.01 


-2.0 ±0.2 




Y 


G 


2131—021* 


ct 


4 


8 


105 


1 7 4- 0'^ 


7.7 ±1.6 


1QQ1 Q4- 1 4. 


N 


Q 


z i J4+UU4 


a ' 


O 


1 o 


— yy 


n n7 4- n n7 
U.Uz it U.Uz 


i.J it 1. / 




I 


1 ■ 

c. 


2136+141* 


B ' 


6 


0.3 


-72 


0.02 ±0.01 


1.8 ±1.4 




N 


G 


2144+092 


Bt 


3 


0.6 


83 


0.03 ± 0.04 


1.5 ±2.0 




N 


F 


2145+067* 


B 


8 


0.3 


-95 


-0.01 ±0.01 


-0.6 ±0.4 




N 


G 




Ct 


8 


0.6 


-119 


0.027 ± 0.003 


1.4±0.2 


1974.2 ±2.3 


N 


G 




D 


8 


1.0 


-118 


0.03 ±0.01 


1.4 ±0.4 


1960.4 ±9.6 


N 


G 


2200+420* 


Bt 


7 


2.4 


-172 


1.41±0.13 


6.5 ±0.6 


1995.1 ±0.2 


Y 


E 




C 


5 


2.5 


-171 


1.12 ±0.22 


5.2±1.0 


1995.6 ±0.5 


Y 


E 




D 


4 


2.7 


-172 


0.99 ±0.18 


4.5 ±0.8 


1996.7 ±0.5 


Y 


E 




E 


3 


1.9 


-161 


1.09 ±0.37 


5.0±1.7 




Y 


P 
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TABLE 2 — Continued 



Object 


Comp. 


N 


<«> 


<i9> 


Mr 






Distinct component 


Rating 








(mas) 


(deg) 


(mas/yr) 




(years) 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 




13 ' 


o 


1 c 
z.K 


— 14j 


U. j4 It U.UZ 


O.J it U.J 


lyoy. 1 it U.O 


I 






nt 


4 


1 c 
1 .J 


1 Q 


U.U4 ±: U.U4 


T 1 J- T /I 

Z.J it Z.4 




IN 




2223-052* 




6 


1.9 


102 


0.49 ±0.09 


32.5 di 6.0 


1994.0 ±0.8 


Y 


E 




C 


7 


5.2 


98 


0.31 ±0.03 


20.4 ±1.9 


1981. 4± 1.6 


Y 


E 


2230+1 14 


B ' 


8 


1.9 


143 


0.03 ± 0.04 


1.4 ±2.1 




Y 


E 






o 


4.0 


1 

loo 


— U.Uj it U.UJ 


T "7 _L 1 A 

— Z. / it i.O 




I 


Ji 




D 


5 


6.5 


156 


0.00 ±0.01 


0.0 ±0.8 




Y 


E 




E 


3 


17.8 


141 


-0.23 ±0.37 


-12.3 ±19.9 




N 


P 


2234+282 


Bt 


5 


1.2 


-97 


0.12±0.05 


5.1±2.2 




Y 


E 


2243-123* 


Bt 


4 


2.7 


15 


0.29 ±0.03 


10.7±1.1 


1988.7 ±0.9 


Y 


E 




C 


3 


10.3 


25 


-0.14 ±0.05 


-4.9 ±1.9 




Y 


G 




D 


4 


11.3 


33 


0.11±0.10 


4.1±3.5 




Y 


G 


2251+158* 


Bt 


8 


5.2 


-86 


0.04 ±0.04 


1.9±2.1 




Y 


E 


2345-167* 


Bt 


5 


1.7 


124 


0.03 ±0.05 


1.0±1.6 




Y 


E 



Note. — Columns are as follows: (1) lAU Name; *Source is a member of the representative sample; ^Motions for 0316+413 which are taken from unpublished data do not appear 
in Figure 1. The motions for 0238—084 are from Vermeulen et al. i 2003a). (2) Component Identifier; ^Component is the brightest. (3) Number of epochs with measured position; (4) 
Mean radial position, relative to the core; (5) Mean structural position angle, relative to the core; (6) Angular radial speed and Icr uncertainty; (7) Radial speed in units of the speed 
of light and Icr uncertainty, computed assuming H{) = 70 kms~' Mpc~', Sim = 0.3, and SIa = 0.7; (8) Extrapolated epoch of origin and uncertainty; (9) Does the component satisfy 
criterion "b" in Section 4.2, i.e., is it a well defined feature?; (10) Overall Quality Rating. 
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TABLE 3 
Vector proper motion results. 



Object 


Comp 


N 


<ff > 


<i9> 




4> 


Papp 


\< '&>-(!>] 








(mas) 


(deg) 


(mas/yr) 


(deg) 


(deg) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(V) 


(8) 


(9) 


no^S-i-ili ^ 




5 


1.5 di 0.0 


1 UO . / ZL. \J.O 


1 -1- 09 
u. i u in u.uz 


191 Q _|_ Q 9 
1 Z 1 .0 HZ O .Z 


6.6 ±1.2 


18.1 ± 8.2 






5 


0.9 di 0.0 


-54.2 ±1.6 


A 1 4. -1- n 09 


—52.3 ± 7.3 


0.2 ±0.0 


1.9 ± 7.5 


\J i \j\J U i J 


c 




9 1 _i_n 


—1 12 7 ± 1 2 


A 97 _|_ A A^ 

u.z / ^ u.uo 


—110/^4-7 4 

1 l7.D ^ / .M- 


22 4 ± 2 1 


Q -1- 7 <; 

D.7 ^ / .O 


nHQ-i-91S 


c 


D 




Zi .7 ^ 1 . i 


A 90 _|_ A Avl 

U.Z7 ^ u.ut 


_9/; 1 _|_ Q Q 

ZD.O ^ 7.7 


1 Q 7 -1-9 1 

lO. / ^ Z.O 


4 C -1- Q Q 
t.O ^ 7.7 


0'?'^ -1- 78^ 




fi 
O 


'I c _|_0 1 


— 1 9 rf; -1- rfi 


9^4-0 0^ 


11 6±4 2 


14 7 ± 2 8 


24 2 ± 4 3 




g 


9 


2.7 di 0.0 


izo.o 31 u.z 


18 ±0 01 

U. lO HI U.Ul 


1 90 1 4- 1 8 

IZU.O HZ 1 .0 


11.0 ±0.4 


3.0± 1.8 




Q 


9 


4.9 di 0.0 


i9(C 0±0 S 

IZU.U ^ U.J 


A 90 _|_ A 09 
u.zu ^ u.uz 


1 o o.u ziz u.o 


12.2 ± 1.4 


9.0 ± 6.8 


O'^^fi — OIQ 


g 


5 


1.4d=0.1 


66.9 ±1.8 


A 99 _|_ A 04. 

U.ZZ HI U.UH- 


67.1 ±7.0 


10.1 ± 1.7 


0.2 ±7.2 


\JH1.J -T- J /y 


J) 


o 


O.J !!n U. 1 


00.7 ^ u.o 


1 S9 4- 04. 

1 .OZ ^ U.UM- 


1 -I- 1 4. 

DO.O 3Z l.M- 


4. Q -1-0 1 

M-.7 HZ U. 1 


ft 4- 1 ft 

U.D ^ 1 .D 


fU'^n 4-0^9 


g 


10 


4.7 di 0.1 


— 1 10 8 ± n 

1 lu.o m U.U 


1 77 4- 06 
1 . / / ^ u.uu 


—108.7 ±1.1 


3.9 ±0.1 


2.1 ± 1.3 




c 


c 
o 


D. 1 ^ U.Z 


— 10*; Q-l- 1 1 

1U0.!7 ^1.1 


1 Q1 _|_o 17 

1 .Oi ^ U. 1 / 


—1 00 4- '^ 7 
lUU.U ^ o. / 


4 04-0 4 
t.u ^ u.t 


0.7 ^ o.o 






Q 

y 


-1- 1 
O.U IC V^, 1 


1 1 9 n -k 1 Q 

— 1 IZ.U It 1 .V 


1 4- A 1 

1 ,OV III U. 1 u 


1 0'^ Q 4- 9 '^ 
1U0.7 ^ Z.O 


1^4-09 
O.o zt U.Z 


^; 1 -1- 1 9 

D. 1 HZ O.Z 






J 




— 1 1 o .+ m z.u 


Z..\)0 HZ U.ZU 


_i 1 4 C _|_ -5 < 
— i IH". O ^ .o 


4 ft 4-0 4 

t.D ^ U.'r 


3 6 ± 4 4 


07'^^ + 1 78 

U / J J X / o 


D 


5 


4.0di 0.1 


63.9 ± 0.9 


IS 4- 0^ 


—119 1 ±126 

J.X7.X ZZ 




1 77 -1- 1 9 7 

1. / 1 .\f HZ / 


07'*fi4-017 

U / T U J. / 


g 


5 


4.2 d: 0.1 


-70.0 ± 1.6 


1/.7IJ _]_ V.lfT' 


—87.3 ± 4.3 


11.9 ± 0.5 


17.4 ± 4.6 








10 di 1 


O / .0 ^ U.Z 


A aq -J_ A AA 

U,07 HI u.uo 


—86 1 ± 3 1 


4 8 ± 8 


1 7 -1- -5 1 
i . / ^ 0. i 






J 


9 1 -1-0 1 


Q _L 9 9 
O / . o III Z . Z 


A Q/; _|_ A AO 


79 -1- 1 8 
— 1 L.yj ^ . o 


1 ft 4- 1 1 
iU.D ^ i . 1 


4 9-1-44 


A7'^S -1-'^ 1 

U / JO J U 


g 


g 


9 8-1-00 

Z.O ^ \J.\) 


1 7q A _|_ n Q 
i /i'.+ m u.o 


A 07 _|_ A 1 
U.U / HI U.Ui 


1 <7 7 _|_ 7 ^ 
10 / . / ^ / .D 


9 S 4- 4 

Z.O ^ U.M- 


91 ft -1- 7 ft 

Z i .O ^ / .D 




Q 




3.1 di 0.1 


—63 1 ± 1 


A J_ A OA 


CO 0± 6 S 
ou.u ^ u.o 


8.2 ±1.0 


13.1 ± 6.5 




g 




1 7 di 


1 20 1 ± 4 


A 19 4-001 

1.7. i Z _I_ 1.7.1.7 1 


129.3 ± 4.5 


5 9 ± 5 


9.3 ±4.5 


07^4 -1- 1 0ft 


g 


5 


4 Odi 1 


14 6 ± 9 


7*; 4- 1 ^ 

u, / o m u. lo 


— 7 1 ± 4 5 


12 9 ± 2 1 


21 7 ± 4 6 


OR 1 4. -1- 4.9 S 


g 


D 


Q -1- f > 


8Q ^1 -1- 8 

o^.u ^ u.o 


18 4-0 09 

U. io ^ U.UZ 


8^1 S 4- 9 
OU.O ^ z.u 


9 A -L 1 
Z,7 ^ U.O 


114-91 

0.1 ^ Z. 1 


n89'^-i-n'^'^ 


c 


O 


9 4. -1-0 1 


14. 1 -1-0 

It.l ^ U.D 


4.8 4- 04. 


10 8 4-1 f^ 
lU.O ^ l.D 


14 4- 1 1 

If .O ^1.1 


O.O ^ 1 . / 




J) 


o 


1 1 -1-0 1 

1.1 !E U. 1 


1.7 ^ O.U 


'11 4-0 0*^ 

U.Ol HI U.UO 


Ol .O ^ D-t 


4.4- 1 S 

7.f ^ l.O 


ft 4-7 1 

U.D HI / . 1 


nsi'^i 4-909 

Wo J 1 T IXjZi 


n 


J 


10-1-0 

l.U !E U.U 


—01 0-1-10 

— 71.7 HI l.U 


^9 4-0 04. 

U.OZ HI U.U^ 


—00 ft 4- 9 ^ 

— 7U.D HI Z.O 


4-0 

7.7 ^ U.7 


19-1-97 
1 .z ^ z. / 




J) 




n Q -1-0 1 


—1 09 4. 4- 1 

— lUZ.M-^ 1.0 


'18 4- Ofi 

U.Oo 31 U.UO 


—111 4 4- 

— 1 1 l.t HI o.u 


7 9 4-1 1 

/ .Z ^ 1.1 


Q 04- 1 

7.U ^ O.O 


0Q17-4-44Q 


g 




1.2 di 0.0 


177.7 ±0.5 


07 ± 01 

U.U / ^ U.Ul 


—176.1 ± 5.7 


5.9 ±0.8 


6.2 ± 5.7 


xyy^o A J** 


Q 


5 


l.OdzO.O 


— 129.3 ± 1.4 


10 ±0 no 

u. J.U ^ u.uu 


—96.6 ± 9.5 


4.2 ±0.1 


32.7 ±9.6 


1019 -i-9'?9 


g 


5 


9 9 _|_ A n 


1 AO 7 _|_ A /I 

iV/7. / ^ V/,f 


97 4-0 09 


11 1 7 ± 1 8 


Q 4- ft 

7.U HI U.O 


9 04-18 
Z . 1. J ^ 1 . o 




g 




1 3 di 1 


—172 3 ± 4 


91 4- 04 
U.ZX It u.ut 


1 77 9 _i_ ft 

A / / .Z ^ U.O 


19 ft 4- 9 9 

AZ.O It Z.Z 


10 4 ± 7 


1101+ 


Q 


7 


1.4 di 0.0 


K ± 1 1 


21 ± 02 


—57.1 ±3.8 


0.5 ± 0.0 


18.3 ± 4.0 




J) 




0.8 di 


an a _L 9 o 
-J (./ . -J HI z . o 


1 7 4- 09 

yj. 1 / HZ 1./.V7Z 


—26.8 ±92 


4 ± 0.0 


3.5 ±97 


1127 — 145 


Q 


7 


4.1 di 0.0 




n 09 + n 01 


70.4 ± 3.1 


5.5 ± 0.8 


12.6 ±3.1 








11 di 1 


1 0/i S -1- 
luu.o m u.o 


«Q -1- 01 
\f.jy HZ u.uo 


141 ft 4- 1 

AH'O.O HZ O.U 


4 ± 2 


Ift 7 -1- 1 

oo. / It o.u 




g 




6.1 ± 0.1 


-112.7 ± 0.4 


1 OS + 04 


-119.6 ± 1.6 


10.9 ± 0.4 


6.9 ± 1.7 






7 


Q A _i_ n 1 

y .o ni u. i 


_i 9fi 9 _|_ n 7 
— izu.z HZ u. / 


1 Ifi -1- 08 
1 ,OU HZ U.V^O 


— 119 1 ± 2 9 


14 1 4- Q 

IM-. i ^ u.7 


1.1 ^ o.u 




iJ 


c 
o 


10 8 -1- f > 9 


— 1 IV, J m 1 .u 


U.o J ZlZ U. iO 


1 1 7 -I- 1 1 

— 1 iV. / ZL iU.O 


8 ft 4- 1 7 
o.u It i . / 


9-1-101 
U.Z It iU.O 






10 


3 3 di 1 


—1 16 1 ± 1 1 


1 97 _|_ 07 

L .iL / in U.U / 


_i 1 c 4 _|_ 9 a 
i 10.+ HZ Z.O 


13 1 ± 8 


7-1-9 S 
V^. / HZ Z.O 




P 




18.7 diO. 2 


_i9a 1 -\-C\ 

iZO . i ^ U.U 


8^ 4- M 

u.oo ^ U. IH- 


_140 4- 10 ft 

IH-U.O ^ lU.U 


8.5 ±1.5 


17.4 ± 10.6 


12S'? — O^S 




10 


3.3 di 0.0 




29 ± A 01 


—1 10 1 ± 1 1 

— 1.JV.J ^ J. A 


9.1 ±0.4 


14.4 ±3.1 


l'^4.S + 19*^ 

1 Of J 1 Z,J 




J 


9 0-1-0 


1^y4. 9 -1-0 8 
luf .z It u.o 


lfi-l-0 01 

U. ID HI U.Ul 


1 7-2 a _L 7 <i 

1 / J.J 1 .J 


19-1-0 1 

l.Z ^ U. 1 


Q 1 4.7 1 
7. 1 HI / .0 






e 


Q S -1- 


11^1 S -1- 1 
XDX.O IC U.A 


04 4- 01 
U.ut It U.Ul 


114 14-15 1 
AA4.J HI A3. A 


14-01 
U.J It U.A 


47 ? 4- 1 5 1 
1/ .3 It A3. A 


1 J lu yjoy 


g 


5 


2.8 di 0.1 


—26.5 ± 0.9 


8S ±0 04 


— 23.2 ± 2.1 


19.0 ±0.8 


3.3 ± 2.3 


1 5'*9 + 01 




/: 
o 


1 1 -1-0 n 


1 4 4- 1 

AJU.t HI A.U 


91 4-0 01 

\J.A± It U.UA 


1 4ft 4- 9 5 


14 14-0 7 
Xt.J ^ u. / 


16 4 ± 2 7 




g 


5 


4.4 di 0.0 


10.4 ±0.3 


1 2 ± 01 

U.Aai ^ U.UA 


74.4 ± 2.9 


7.7 ±0.8 


64.1 ±2.9 


1 ^lO^i + 1 0^1 

i DUD i UO 


g 


5 


9 dr 


CO A _L 9 C 
— J o ,U HI Z.O 


Oil 4- 09 

yj.J 1 It V/.UZ 


_47 Q _|_ 4 
/ .o It '-*.\) 


10 7-1-1 1 

i O. / It 1 .o 


10 2 ±4 7 


1 A'i'i + '?S9 

1 DO J JOZ, 


g 




1 S -1- 


_g A 9 _|_ A S 
— oU.Z HZ U.O 


1 S 4- 01 
u, io It u.uo 


_8ft 9 4-10 
oD.Z HI Ji.y) 


11 4 4-1 Q 

i 1 .H- It i .7 


() ± 3 1 


1 Ut i T J ^ 7 


g 


9 


1.7 di 0.0 


70 .U ^ V^.U 


4Q 4- 09 


—99 9 ± 1 7 


1 7 _[_ § 


5.0 ±1.8 








Q -1-0 
yj.y ^ U.U 


_1 AA A _L A 7 
iUU.D 1 


17 4- 01 
U.O / HZ U.UO 


—92 9 ± 2 5 


1 9 q -L 

i Z , 7 It \j.y 


7 7 ± 2 6 


1 ^i4.9 + ^^QO 


g 




2 2 dr 1 


— IDZ.H-HZ i.U 


n 4-0 04 

U.Oo HZ U.UH- 


—1 ft 1 4-48 

i O i . U It , o 


1 6 ± 1 6 


1 4 ± 5 


l^;'^S+077 


g 


J 


3 7 di 


— Oo,+ HZ U.U 


41 4- 01 

U.H-O HZ U.UO 


—4ft ft 4- 1 4 

tU.O H! O.H- 


15 4 ± 1 


8 2 ± 3 5 






5 


k 04- 
o.u in U.U 


—44 2 ± 2 


1 4 -1- 09 

U.At It U.UZ 


17 1 ± 6 1 


4 4-00 

■+.7 HZ U.7 


ftl 1 4- ft 1 

DA.O Zt D, A 


1 SOO + 440 

i OUU \ tH-Vj 


g 


5 


2.6 di 1 


— 163 1 ± 1 4 


40 4- 04 
yj.^y HZ u.v^T- 


—155 7 ±4 3 


18 6 ± 1 6 


7 4 ±4 5 


1 84^ + 7Q7 


g 




5 5 di 


—40 9 ± 3 


SS 4- 01 

\t.JJ HZ U.UO 


_9Q 0-1-9 9 

Z7.7 HZ Z.Z 


9 1 -k 1 

Z, A It U.A 


10 9 ±2 2 




c 


10 


7.2 di 0.0 


± 2 


61 + 02 


-26 6 ± 1 


2.3 ± 0.1 


9.3 ± 1.0 




g 


c 
o 


9 1 -1-0 1 

Z. i ^ U. i 


^-1-0 7 
iDj. J in U- / 


90 4- 01 

U.Z7 ^ u.uo 




S ft -U ft 
O.D ^ U.D 


7 04-1 1 

/ .U ^ 0.0 


1Q';7+4.0S 


J) 


e 
o 


■^0-1-0 
o.u ^ U.U 


/o.D in u.^ 


1 8 4- 09 
u. lo ^ u.uz 


_70 9 -L 9 n 


7 4-0 1 

U. / ^ U. 1 


44-9 
u.f ^ z.u 






g 


2.2 di 0.0 


—79.7 it 0.4 


1 Q 4- 03 

U. 17 ^ U.UO 


— 75.6±3.1 


0.7 ±0.1 


4.1 ± 3.1 




F 


6 


1.7ib0.1 


-79.4 ±0.5 


0.22 ±0.03 


-79.5 ±1.9 


0.8±0.1 


0.1 ±2.0 




G 


8 


1.2=t0.1 


-79.0 ±0.7 


0.16 ±0.03 


-80.2 ±2.9 


0.6 ±0.1 


1.2 ±2.9 


2005+403 


B 


5 


3.2 ±0.1 


119.8 ±1.4 


0.20 ±0.03 


83.6 ±13.6 


14.9 ±2.1 


36.2 ±13.7 


2113+293 


B 


5 


0.3 lb 0.0 


177.0 ±2.2 


0.04 ±0.01 


-131.3 ±8.2 


2.4 ±0.3 


51.7 ±8.5 


2145 + 067 


C 


8 


0.6 lb 0.0 


-119.4 ±0.7 


0.03 ±0.01 


-156.4 ±5.3 


1.8 ±0.3 


37.0 ±5.3 


2200+420 


B 


7 


2.4=b0.1 


-171.6±1.6 


1.41 ±0.12 


-170.7 ±3.1 


6.5 ±0.6 


0.9 ±3.5 




C 


5 


2.5 ±0.2 


-171.4 ±1.5 


1.16 ±0.22 


172.8 ±3.4 


5 3±1.0 


15.8±3.7 


2201+315 


B 


6 


2.8 ±0.0 


-144.5 ±1.2 


0.34 ±0.03 


-142.7 ±5.3 


6.3 ±0.5 


1.8±5.4 


2223-052 


B 


6 


1.9±0.1 


102.4 ±2.8 


0.50 ±0.08 


107.3 ±6.4 


32.9 ±5.4 


4.9 ±7.0 




C 


7 


5.2 ±0.1 


98.0 ±0.3 


0.31 ±0.03 


90.6 ±3.3 


20.6 ±1.8 


7.4 ±3.3 



Note. — Columns are as follows: (1) lAU Name (B1950); (2) Component Identifier; (3) Number of epochs with measured position; (4) Mean radial position, relative to the core; 
(5) Mean structural position angle, relative to the core; (6) Angular velocity; (7) Direction of velocity; (8) Linear speed in units of the speed of light; (9) Misalignment angle. Sources 
with significant non radial motion are shown in bold print. 



